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This paper presents adaptive boundary element methods for positive, negative, as well as zero order 
operator equations, together with proofs that they converge at certain rates. The convergence 
CO rates are quasi-optimal in a certain sense under mild assumptions that are analogous to what is 

' typically assumed in the theory of adaptive finite element methods. In particular, no saturation- 

,——1 type assumption is used. The main ingredients of the proof that constitute new findings are some 

■^^ results on a posteriori error estimates for boundary element methods, and an inverse-type inequality 

►^ involving boundary integral operators on locally refined finite element spaces. 

g 1 Introduction 

Let r be a (closed or open) polyhedral surface in M^ . We consider equations of the form 
^ Au = f, (1) 

^^ where f £ H *(r), and A : H^{r) — )• H *(r) is an invertible linear operator. Strictly 

speaking, the function spaces should be slightly modified if F is an open surface or t > 0, but 
for the sake of this introduction we will gloss over this point. The operators of interest are 
the boundary integral operators that arise from reformulations of boundary value problems 

^—^ as integral equations on the domain boundary, cf. McLean (2000). Then the problem (1) 

J> corresponds to a boundary integral equation, and a very popular class of methods for its 

'k> numerical solution is boundary element methods (BEM), which can crudely be described 

^ as finite element methods applied to boundary integral equations, cf. Sauter and Schwab 

^ (2011). As with finite elements, there is the adaptive version of BEM, whose main feature 

is to automatically distribute mesh points hopefully in an optimal way so as to obtain an 
accurate numerical solution. Even though those methods perform very well in practice, 
their mathematical theory is in a not very satisfactory state, especially if one compares it 
with the corresponding theory of adaptive finite element methods. In the latter context, the 
sequence of papers Dorfler (1996), Morin, Nochetto, and Siebert (2002), Binev, Dahmen, 
and DeVore (2004), Stevenson (2007), and Gascon, Kreuzer, Nochetto, and Siebert (2008) 
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has laid a fairly satisfactory foundation to mathematical understanding of adaptive finite 
element methods for linear elliptic boundary value problems. Specifically, it was established 
that standard adaptive finite element methods generate a sequence of solutions whose error 
decreases geometrically, and that the number of triangles in the mesh grows with an optimal 
rate. Moreover, rigorous treatments of numerical integration and linear algebra solvers are 
within reach. 

However, there has been a significant gap in the current mathematical understanding of 
adaptive BEM. The first steps toward closing the gap have been taken in Carstensen and 
Praetorius (2009), where convergence is guaranteed for an adaptive BEM with a feedback 
control that occasionally adds uniform refinements, in Ferraz-Leite, Ortner, and Praetorius 
(2010) where geometric error reduction is proven under a saturation assumption, and in 
Aurada, Ferraz-Leite, and Praetorius (2011), where convergence is established under a 
weak saturation-type assumption. To give an idea of what a saturation assumption is, 
one form of it requires that if Uk is the numerical solution of (1) at the current stage of 
the algorithm execution, and if Uk is the would-be numerical solution had we replaced the 
current mesh by its uniform refinement, then \\u — UkWH^/r) ^ (3\\uk — Uk\\H^(r), where 
/3 > is a constant. Such assumptions were common in finite element literature before 
Dorfler (1996) and Morin et al. (2002) proved geometric error reduction without relying 
on a saturation assumption. 

In this paper, we prove geometric error reduction for three kinds of adaptive boundary 
element methods for positive (t > 0), negative (t < 0), as well as zero order operator 
equations, without using a saturation-type assumption. In fact, several types of saturation 
assumptions follow from our work as a corollary. Moreover, bounds on the convergence 
rates are obtained that are in a certain sense optimal. 

During the final stage of this work, and independently of it, a preprint Feischl, Karkulik, 
Melenk, and Praetorius (2011) appeared in which the authors announced a proof of certain 
optimality result for an adaptive BEM for a negative order operator equation, based on the 
a posteriori error estimator from Carstensen, Maischak, and Stephan (2001). Details are 
not given, but the general outline shows some similarity with our main line of reasoning. 

A large chunk of the techniques developed in the adaptive finite element theory is not 
specific to differential equations, providing a nice starting point for us and a clue on what 
ingredients are missing in the boundary element theory. In this work, we needed to supply 
two kinds of ingredients that constitute new findings: some results on a posteriori error 
estimators for BEM and an inverse-type inequality involving boundary integral operators 
and locally refined meshes. The issue with the theory of error estimators has been most 
obvious; what is usually guaranteed is only one of the two inequalities that are necessary 
for a quantity to resemble the error. There is a very few estimators with both upper and 
lower bounds proven; one can mention the estimators proposed in Faermann (1998, 2000, 
2002), see Carstensen and Faermann (2001) for a more thorough discussion. To name 
some of the relatively recent works on this subject, we have Carstensen, Maischak, and 
Stephan (2001), Carstensen, Maischak, Praetorius, and Stephan (2004), and Nochetto, von 
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Petersdorff, and Zhang (2010), where upper bounds for certain residual- type estimators are 
estabhshed, and Erath, Ferraz-Leite, Funken, and Praetorius (2009a), and Erath, Funken, 
Goldenits, and Praetorius (2009b), where upper and lower bounds are proven for a large 
number of non-residual type estimators, with the upper bounds depending on various forms 
of saturation assumptions. Now, even if one had both upper and lower bounds for an error 
estimator, one still needs a so-called local discrete bound before attempting to apply the 
general techniques from the finite element theory. Such an estimate has been entirely open 
for boundary element methods, except the recent claim in Feischl et al. (2011) regarding 
the residual-type estimator from Carstensen et al. (2001). 

In this work, we establish all missing bounds for a number of residual-type error esti- 
mators for positive, negative, as well zero order boundary integral equations, including the 
estimators from Carstensen et al. (2001), Carstensen et al. (2004), and Faermann (2000, 
2002). Some of our bounds contain the so-called oscillation terms, that give useful esti- 
mates on how far the current mesh is from saturation. Note that analogous terms also arise 
in the finite element theory. In order to handle the oscillation terms, which turned out 
to be much more involved than the corresponding affairs in the finite element context, we 
prove an inverse-type inequality involving boundary integral operators and locally refined 
meshes. A special case of this inequality has been announced in Feischl et al. (2011). Our 
proof of the inverse- type inequality in general requires the underlying surface F to be C^'^ 
or smoother, but for open surfaces it allows the boundary of F to be Lipschitz. So in 
general, polyhedral surfaces are ruled out, which is likely an artifact of the proof, and will 
hopefully be rectified in the near future. 

This paper is organized as follows. In the next section, we fix the general setting of 
the paper, and recall some basic results that will be used throughout the paper. Then 
in the following three sections, namely in §3, §4, and §5, we consider adaptive BEMs for 
zero, positive, and negative order operator equations, respectively. In each of the three 
cases, we first study certain residual based a posteriori error estimators, then design an 
adaptive BEM based on those estimators, and finally address the question of convergence 
rate. There are some general arguments and remarks that can be applied to all of the three 
cases, and so in order not to be too repetitive, they shall be discussed in §3 for the zero 
order case in detail, and then in §4 and §5, we will simply refer to them if needed. The entire 
analysis depends on an inverse-type inequality involving boundary integral operators and 
locally refined meshes, which then is verified in §6 for a general class of boundary integral 
operators. We will conclude the paper by summarizing the results and making a list of 
important open problems. 

2 Generalities 

In this section, we will set up the necessary vocabulary and collect some basic results that 
will be used throughout the paper. Let fi be a compact closed, n-dimensional, patchwise 
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smooth, globally C^^^'^ manifold, which will be the habitat of all functions and distribu- 
tions that ever occur in this paper. In practice, we typically have n = 1 or n = 2, so in the 
discussions that follow we often use the language of n = 2, and we shall indicate whenever 
we lose generality by implicitly restricting to the case n = 2. We will assume that each 
smooth (closed) patch of fi is diffeomorphic to a polygon, so one can think of fi as the 
surface of a bounded polyhedron, with faces and edges now allowed to be smooth surfaces 
and curves, respectively. For s G [— z^, z^], we denote by H'^{il.) the usual Sobolev space of 
order s on fJ. Let a; C f2 be an open subset of ri with Lipschitz boundary. Then we define 
the following two kinds of Sobolev spaces 

H'{lo) = {u\^ : u G H'iQ)}, and H'iuj) = {u e H'iQ) : suppw C UJ}, (2) 

for s e [-i^,i^]- Obviously H%Q) = H^'^Q), and it is known that H%uj)* = H~''{uj), 
cf. McLean (2000). Note also that H^{u:) is a closed subspace of H^{uj) for s > 0. The 
definitions (2) give rise to the canonical norms on H^{uj) and H^(uj) inherited from the norm 
on H^{Q), and at least for s > 0, these norms are known to be equivalent to certain norms 
defined either by interpolation, or in terms of moduli of smoothness, with the equivalence 
constants depending only on the dimension n, the order s, the Lipschitz constant of cj, 
and the particulars (i.e., the local coordinate patches and the partition of unity) in the 
definition of H^{i},). Let us expand on this a bit. Since the spaces with s > 1 are of 
secondary importance to us, for simplicity here we focus only on the case s £ [0, 1], hence 
by duality on \s\ < 1. With || • ||^ and | • li^^j denoting the L^-norm and the usual H^- 
seminorm on uj, respectively, we define | • iH^iu}) to be the interpolatory seminorm between 
II • \\uj and I • |i^^, for concreteness by using the iC-functional. Then || • ||/^s(^) := | • \h={u)) is a 
norm on H'^{u}), which can be made into a norm || • ||j|/s(j^-) on H^{u}) by combining it with 
the L^-norm. It can be shown that these norms are equivalent to the canonical norms on 
H^{u}) and respectively H^{uj), with equivalence constants depending only on n, s, and the 
particulars of the definition of H^{Q), cf. McLean (2000). Moreover, we have the following 
useful property that if wi, . . . ,0;^ are disjoint Lipschitz domains such that IJi'^« — ^-^ then 
for |s| < 1 we have 

i i 

for u G H^{(jj) and v G H^{uj), with the constant Cg > depending only on s, and the 
norms for s < defined by duality, cf. Ainsworth, McLean, and Tran (1999) and von 
Petersdorff (1989). Another interesting norm on H^{u}) for s G (0,1), is the Slobodeckij 
norm 

||„,||2 _||„,||2 , I |2 -^1 I |2 _ f f k(^) -Viy)\ . . (A^s 
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where \\v\\uj denotes the L^-norm on w. It is immediate that if wi, . . . ,Ljfc are disjoint and 
Lebesgue measurable sets such that |Jj uji C uj, then for s > we have 



El 



"l.„ < l"ll:.. ^ e -W'H. (5) 



which is the counterpart of the first inequality in (3). The seminorms | • \h'(ui) ^-nd | • |s,a; are 
equivalent, with the equivalence constants depending only on n, s, the Lipschitz constant 
of u, and as usual the particulars of the definition of H^{i}). This fact can be proven as 
in McLean (2000) by relating the both seminorms to the canonical norms defined through 
the norm on H^{Q). A more direct way to relate the two seminorms would be to first 
connect the Slobodeckij seminorm with the Besov-style seminorm defined by moduli of 
continuity, and then use the equivalence between the moduli of continuity and the K- 
functional, established in Johnen and Scherer (1977). Note that one has to inspect (and 
adapt) the proofs in McLean (2000) and Johnen and Scherer (1977) to reveal the relevant 
information on the equivalence constants. It should be emphasized that since we shall 
be dealing with an infinite collection of domains, in partiuclar of sizes that can shrink to 
zero, when using norm equivalences one must be careful about ensuring a control over the 
equivalence constants. In the current setting, the equivalence constants do not depend on 
the size of uj, and the Lipschitz constants are controlled by restricting the class of domains 
to the so-called shape regular triangles, see below (14). In the proofs, we make an effort to 
use the interpolatory norms for as long as possible, and so to apply the norm equivalence 
only when necessary. 

In the following, we fix L C $7 to be the whole of il, or a connected open set whose 
boundary consists of curved polygons. This will be the domain on which we consider our 
main operator equation 

Au = f, (6) 

where / S H~^(T), and A : H^{T) — )• H~^{T) is a linear homeomorphism. We assume that 
the operator A is self-adjoint and satisfies 



{Av,v)>a\\v\\^,^^y \\Av\\H-t(r) < PWvWfjt^Y'), (7) 

for V G i?*(r), with some constants a > and /? > 0, where (•, •) is the duality pairing 
between H~^(r) and H^{r). We introduce the energy norm ||| • ||| = {A-, •)^", and note that 
it is equivalent to the i7*(r)-norm 



/all • ll^.(r) < III • III < V/3II • htiry («) 

We also have the norm equivalence 

all • llj^t(r) < P • llH-*(r) < /5|| • lli?t(r)' (9) 

which is the basis of all residual based error estimation techniques. 
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Suppose that a closed linear subspace S C -ff*(r) is given. Then the Galerkin approxi- 
mation Us ^ S of u from the space S is characterized by 

{Aus,v) = {f,v), yveS. (10) 

We have the Galerkin orthogonality 

III III 2 I III III 2 III III 2 vj ,^ ct /'i -1 ^ 

|n — U5I + |ii5 — ul =|n — 1;| , vv € b, (llj 

which implies that us is the best approximation of u from S in the energy norm, and that 
the Galerkin approximation is stable: 

V^hsWiit^r) ^ i^slll < |||m||| < -j^UWh-^t)- (12) 

In light of (9), the residual rs = f — Aus is equivalent to the error: 

"Ik - ^sllHt(r) < lks||H-*{r) < P\\u - usWjjtf^r)- (13) 

In the sense that the residual is a computable quantity that gives bounds on the true error 
in terms of this equivalence, the first inequality in (13) is an example of a global upper 
bound, while the second one is that of a global lower bound. Upper and lower bounds in 
this context are also called reliability and efficiency, respectively. 

We study the Galerkin approximation by piecewise constant or continuous piecewise 
linear functions on adaptively generated triangulations of the manifold T. Let us now fix 
some notations and terminologies related to this discretization. An open subset of T is 
called a (surface) triangle if its closure is diffeomorphic to a flat triangle, and the latter is 
said to be the reference of the former. Assuming that a choice is made of a reference for 
each surface triangle, notions related to flat triangles can be planted onto surface triangles 
through their references. For instance, in the following, straight lines, midpoint, etc., 
should be understood in terms of the reference triangles. We call a collection P of surface 
triangles a partition of F if F = IJreP ^' ^^'^ t Da = for any two different T,a £ P. For 
reflning the meshes we mainly use the so called newest vertex bisection algorithm, which 
we describe now for the reader's convenience. General discussion on this algorithm can 
be found e.g., in Binev, Dahmen, and DeVore (2004); Stevenson (2008). We assume that 
with any triangle r comes its associated newest vertex v{t), so that when it is needed to 
be reflned, r is subdivided into two triangles by connecting v(t) with the midpoint of the 
edge opposite to it. The midpoint used in the bisection is now the newest vertex of the 
both new triangles. The 2 new triangles so obtained are called the children of r, and the 
reflnement of r is just the collection of its children. The children of a triangle inherit the 
reference map from their parent. A partition P' is called a reflnement of P and denoted 
P ^ P' if P' can be obtained by replacing zero or more r E P by its children, or by 
a recursive application of this procedure. The newest vertex bisection in 1 dimension is 
straightforward, and an extension to higher dimensions is treated in Stevenson (2008). 
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An alternative to the newest vertex bisection is the red refinement procedure, which 
subdivides a triangle into 4 smaller triangles, by connecting the midpoints of the edges 
with each other. In general, a refinement procedure is a map that assigns to a triangle a 
collection of disjoint triangles whose union is the original triangle (up to taking closures). 
The newly created triangles inherit the reference map from their parent. We will assume 
that there is a universal constant Aq < 1 such that the maximum diameter of the children 
of a triangle is at most Aq times the diameter of the parent. This assumption is easily 
satisfied for the newest vertex bisection and red refinement procedures. Although one can 
imagine many different kinds of exotic refinement procedures, we are going to pretend that 
we have only the afore- mentioned refinement procedures, and their finite powers (e.g., one 
can perform the red refinement twice on a triangle and call it one refinement). The red 
refinement is algorithmically the bit simpler one, and a rule of thumb is that we can use 
the red refinement if we may allow discontinuous functions in our approximating subspace 
S, and we need the newest vertex bisection if we want only continuous functions in S. 
However, see Stevenson (2005) for an approach to use the red refinement procedure in the 
context of continuous approximating subspaces. 

A partition P is said to be conforming if any vertex f of a triangle in P is a vertex 
of all r E P whose closure contains v. Throughout this paper we consider only partitions 
that are refinements of some fixed conforming partition Pq of F. We require that the 
references of the initial partition Pq be so that for any pair of surface triangles that share a 
common edge, the parameterizations from the reference triangles to the common edge are 
equal up to the composition with an affine map. The motivation for this is that we want 
the refinements on both triangles to agree on the common edge. A choice of refinement 
procedure immediately leads to the set [Pq] of all partitions that are refinements of Pq- We 
assume that the family [Pq] is shape regular, meaning that 

a, = sup <^ — f- : T E P, P G [Po] ^ < oo, (14) 

I vol(r) J 

where hr = diam(r). Both the newest vertex bisection and the red refinement procedures 
produce shape regular partitions. The set [Pq] is too large in the sense that often we are 
interested in a certain subset of it that has a good property, e.g., we want to single out the 
conforming partitions from [Pq]. Exactly what subset we want depends on the particular 
setting, and at this level of generality we simply assume that there is a subset adm(P) C [Pq] 
called the family of admissible partitions, which is graded (or locally quasi-uniform, or have 
the K-mesh property), i.e., 

CTg = sup <-^:a,TeP,anT^0,Pe adm(Po) > < oo. (15) 

For example, if n > 2 then the conforming refinements of Pq produced by the newest vertex 
bisection are locally quasi-uniform. If n = 1 or if the red refinement is used, one can define 



2 Generalities 



the admissible partitions to be the ones for which the quantity under the supremum in 
(15) is bounded by a fixed number. Finally, note that the shape regularity and local quasi- 
uniformity together imply local finiteness, meaning that the number of triangles meeting 
at any given point is bounded by a constant that depends only on Us, Cg, and n. 
Given a partition P G [Po], we define the piecewise polynomial space Sf, by 

S^ = {ue L2(r) : u G C{r) if d > 0, and u\r € P^ Vr G P} , (16) 

where P^ denotes the set of polynomials of degree less than or equal to d. Note that 
for curved triangles, polynomials are defined through the reference triangles. We have 
Sp C H''^{T) for s < (i+ 2 with \s\ < u. Oi interest to us are only the spaces Sp of piecewise 
constants and Sp of continuous piecewise linears. We will also employ a slight variation of 
Sp, that is the space Sp of piecewise affine functions that vanish on the boundary of F. 

Now we collect some estimates relating different Sobolev norms for finite element spaces 
and their complements. We will indicate if the constants involved in the estimates depend 
on parameters (such as ds) other than d and n. First of all, we recall Faermann's estimate 

||Hl',r <CfY. l^lLw ^ ^ ^'(r), V ±L2 5|,, {s G [0, 1]), (17) 

zeNp 

for all admissible partitions P G adm(Po)) with Cp = CF{as,crg), cf. Faermann (2000, 
2002). Here Np is the set of vertices in P, and uj{z) = int |J r gp.^g-i r is called the 
star associated to the vertex z, with "int" denoting the interior. The same estimate for 
interpolatory norms has been established in Carstensen, Maischak, and Stephan (2001), 
by a very flexible technique. We will be using their technique on several occasions in §5. 

For a partition P G [Pq], let hp G Sp be such that hp{x) = h-r for x £ t £ P. 
We introduce the space H'^{r,Po) for r > as the space of functions v G L'^{Q,) with 
v\r G H'~{t) for every triangle r G Pq. Now for v G H''{T, Pq) with s G [0, 1], let vp G Sp 
be the L^-orthogonal projection of v onto Sp. Then we have the direct estimate 

\\v-Vp\\r<Cjhl\v\s,r, T G P, (18) 

where the constant Cj = Cj{as)- An immediate consequence is that 

\\hUv-vp)f^<CJY,hl^^+^^\v\l^, (19) 

for Q C P, 7 = int UreQ ^' ^^^ r gR. By using a duality argument, and the bounds (19) 
and (5), one can show also that 

\\v-vp\\ji.s^^^<Cj\\h'pv\\^, veL^T). (20) 

For the continuous piecewise linears, the L-^-projection is nonlocal, and so for conve- 
nience we will employ a quasi- interpolation operator Qp : L^(r) — )• Sp that satisfies 

\\v-Qpv\\r,^ <Cj [maxhr] \v\s,uiM, veH^{T), (21) 

\t6Q / 
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for < r < s < 1 and P G adm(Po)i where oj{'^) = int |J r gp-^^^ /g,! a, and if n > 2, 
adin(Po) is understood to be the conforming partitions created by newest vertex bisections 
from Pq. Recall also that S\, is the subspace of S\, with the homogeneous boundary 
condition. By a quasi-interpolation we mean that {Qpv)\t = (Qp'L'L(t))|t- Examples of 
such operators are constructed, e.g., in Clement (1975); Scott and Zhang (1990); Oswald 
(1994); Bernard! and Girault (1998). Accordingly, in this setting, the estimate (19) is 
replaced by 

\\hUv - Qpv)f^ <C]Y. hl^^'^Ml.iry (22) 

We stated the estimates (18) and (21) in terms of the Slobodeckij norms, but the usual 
way to derive these estimates is by interpolation and norm equivalences, so in particular 
the same estimates hold with interpolatory norms. 
Let us also recall the inverse estimates 

\\v\\Hs[^)<CB\\h~p'v\\^, V £ Sp, (23) 

for s €[0,d+l)n [0, u], and 

\\hpv\\^<CB\\v\\^-s(^^-j, v£SJ,, (24) 

for s € [0, z^], which hold for admissible partitions P € adm(Po)) with the constant Cp = 
CB(o's,o'g). Recall that 7 = intlJ^^QT for some Q ^P. The both inequalities are proved 
in Dahmen, Faermann, Graham, Hackbusch, and Sauter (2004) in a more general setting, 
with a piecewise affine mesh-size function hp. Nevertheless, by local quasi-uniformity, their 
results immediately imply (23) and (24) with piecewise constant /ip, since in our setting 
hp enters only in a derivative-free fashion. Finally, we will make crucial use of inequalities 
of the type 

^ hl^^+^AvWlr < CA\\v\\%^^y V E Sf,, (25) 

that is assumed to hold for admissible partitions P £ adm(Po)5 with Ca = C'a(^) <7s,o"g). 
This inequality is somewhat more demanding than the standard inverse estimates since it 
involves the non-local operator A, and in some sense it requires A to be almost local. In 
fact, boundary integral operators have certain locality property, which is exploited in our 
analysis only through this inequality. Feischl, Karkulik, Melenk, and Praetorius (2011) 
announces a proof of (25) for s = 1, and A equal to the simple layer potential operator 
on a closed surface (hence t = —5)- We prove it in Theorem 6.1 for a general class of 
boundary integral operators. 

3 Operators of order zero 

In this section, we focus on the case where the operator A is of order zero, i.e., the case 
t = 0. This is a nice model case to test our arguments on. On a practical side, this case 
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is a representative of Hilbert-Schmidt operators on general domains or manifolds, and a 
model of boundary integral operators associated to the double layer potential. In the latter 
case, generally one only has Garding's inequality instead of the strict coercivity (7), but we 
expect that such cases can be handled at the expense of requiring a sufficiently fine initial 
mesh, in the spirit of Mekchay and Nochetto (2005) and Gantumur (2008). 

The surface T can be either closed or open. We employ the piecewise constants Sp, 
which we denote simply by Sp. The Galerkin approximation of u from Sp is denoted by 
Up G Sp, and the corresponding residual — by rp. Recall the notations || • || and || • \\ui for 
the L^-norms on F and w C F, respectively. 

The equivalence (13) provides the convenient starting point 

a\\u — up\\ < \\rp\\ < /3||n — up\\, (26) 

which suggests us to use local L^-norms of the residual as error indicators. Below we prove 
a localized version of (26) with the error replaced by the difference between two Galerkin 
approximations. The simple observations in its proof are the essence of this paper, in the 
sense that the rest of the paper can be thought of as an attempt to exploit their natural 
consequences and to extend the arguments to non-zero order operators. Note that no 
explicit condition whatsoever is imposed on the locality of A. 

Lemma 3.1. Let P, P' G [Pq] he partitions with P < P' , and let F* = int UtgpVP' ^' ^^^^ 
we have 

a\\up — upi\\ < ||rp||r* < P\\up — upi\\ + 2\\rp — v\\y*, (27) 

for any function v € Sp' . 

Proof. Recall the definition P\P' = {t G P : t ^ P'}, so F* is the region covered by the 
refined triangles. Let v = upi — up, and let vp G Sp be the L^-orthogonal projection of v 
onto Sp. Then from the Galerkin condition (10) we get 

{Av,v) = {rp,v) = {rp, v — vp) < ||rp||r*||f — fp||r* < lkp||r*||^||r*i (28) 

where we have used that v = vp outside F*. This proves the first inequality in (27). 
To prove the second inequality, let v S 5p/ be supported in F*. Then we have 



l^llr* = ("^j v) = {v — rp, v) + {A{upi — up), v) 
< {\\v - rp||r* + \\A{up' - Mp)||r*) ||^||r* 



(29) 



implying that 



Ikpllr* < Wp — v\\t* + \\v\\t* < 2||rp — v\\y* + ||^(np/ — up)\\, (30) 

which gives the second inequality in (27). □ 
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Remark 3.2. The arguments used in the preceding proof are of course inspired by the 
corresponding finite element theory. However, especially the argument (29) for the lower 
bound seems to have a new flavour, in that it does not break the action of A up into element- 
wise (or star-wise) operations, making it particularly suitable for nonlocal operators. 

We recognize the term ||rp — t;||r* in (27) as an oscillation term, which measures how 
much of the residual is captured when we move from P to P' . It is of interest to control 
this term. Since Sp C H''"(T) for r G (0, |), assuming that A : H^{r) — ?• H^{T) is bounded 
for all r E (0, g), we have Aup G H'"(T) for r in the same range. Now, assuming that 
/ G H'^(r, Pq) for some r G (0, ^), this ensures rp G H^{T, Pq), therefore from (19) we have 

inf \\rp-vfr*<CJ V /if|rp|2^. (31) 

TeP\P' 

This suggests us to define the residual oscillation 

osc{v,P,u):=l Y. hl^\f-Av\lJ , (32) 

\tGP,tCuj J 

for w C r and v G Sp, so that (27) implies 

a\\up — up'W < ||rp||r* < P\\up — up'\\ + 2Cjosc(np, P, F*). (33) 

If the oscillation term is sufficiently small, the difference between two discrete solutions is 
completely controlled by a local L^-norm of the residual. In this sense, the first inequality 
in (33) is an example of a local discrete upper bound, while the second one is that of a 
local discrete lower bound. These bounds are also called local discrete reliability and local 
discrete efficiency, respectively. 

Remark 3.3 (Non-residual estimators and the saturation assumption). Following the wis- 
dom from finite element theory, in this section we will exploit the bounds (33) by designing 
an adaptive algorithm that refines the triangles where the quantity H^pUt- + osc(mp,P, r) 
is large. Before delving into this let us remark that there exist non-residual type estima- 
tors that have some advantage over the residual based estimators. For these estimators, 
usually only the global lower bound is justified, and the global upper bound depends on 
the saturation assumption. The latter is the assumption that fu — up\l < afu — upf for 
some constant a G (0,1), where iip is the Galerkin approximation from some enriched 
space Sp D Sp, which is typically the piecewise constants on a uniform refinement of P, or 
piecewise linears on P. Note that by Galerkin orthogonality, this assumption is equivalent 
to III Up — Up IP > (1 — a'^)|||M — np|p. It is possible to extract some information from Lemma 
3.1 on these issues. For instance, the discrete lower bound in (33), combined with the 
upper bound in (26) gives the second inequality in 

|||up — Mp'l < |||n — up||| < ^|||np — np/||| + 2Cjosc(up,P, F), (34) 
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while the first inequality follows from the Galerkin orthogonality, where ^ > is a constant. 
Taking P' to be the uniform refinement of P, this gives global lower and upper bounds 
for the h — h/2 estimator fup — up'f, and moreover justifies the saturation assumption up 
to an oscillation term. Local bounds for non-residual estimators remain as an important 
open problem. 

In the following, we will fix some r E (0,1], and assumc"'^ that / E H^{T,Pq). Note 
that this is in the same spirit as assuming / E L^ in the context of second order elliptic 
equations, even though there the weak formulation is well-posed for / E H~^. We will 
use the convenient abbreviation osc(P, cj) = osc{up,P,uj). The next lemma collects crucial 
properties of the oscillation osc(P,r) and the combination |||u — up|p-|-osc(P, F)^, assuming 
an inverse-type inequality, which shall be verified in §6. As it turns out, this inverse-type 
inequality is also sufficient to guarantee the finiteness of oscillation (32). The quantity 
|||n — up|p -|- osc(P, r)^, the counterpart of the total error in Gascon, Kreuzer, Nochetto, 
and Siebert (2008), will be the main character in our subsequent analysis. 

Lemma 3.4. Assume that 

^h^^-\Av\l,<CA\lvf, veSp, (35) 

for P E adm(Po); with the constant Ca = C a{A, a s-, (t g) . Then the oscillation (32) is finite 
for any v E Sp and P E adm(Po)- Moreover, the followings hold. 

a) There is a constant Cg > such that 

ju - upf + osc{up, P, Tf < Cg inf (l-u - vf + osc(?;, P, Vf) , (36) 

v^Sp 
for any P E adm(Po)- 

b) There exists a constant A > such that 

osc(P',r)2 < {l + 6)osc{P,rf -X{l + 6)osc{P,r*f + Cs\lup-up'f, (37) 

for any P,P' E adm(Po) with P :< P' , and for any 5 > Q, with Cs depending on 5, 
where T* = int UTeP\P' ^■ 

c) Let P,P' E adm(Po) be such that P :< P' , and let T* = int |J^gp,v pT. If, for some 
/i E (0, 2) it holds that 

\\u-up,f + osc{P',Tf <^i{lu-upf + osc{P,Tf) , (38) 

then with = !^,,\ ~„Jt\ , we have 

\\rp\\l. + osc(P, T*f > 9 {\\rpf + osc(P, r)^) . (39) 



^ In view of (32), one has a computational advantage if r = 1, since there would be no fractional norms. 
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Remark 3.5. The estimate (35) is proved in Section 6 for a general class of singular integral 
operators of order zero, under the hypothesis that A : H^{Vt) — t- H°'{0,) is bounded for 
some a > ^ and that r G (0,o"), where A is an extension of A to $7 if F is open, and 
A = A if r = ^. Recalling that fi is a ^"^'^-manifold, the condition on a translates to 
I' > '^, see Remark 6.3 for details. Therefore Lipschitz curves are allowed, but for n = 2 
we need a C^'^ surface. Note that there is no condition on the regularity of the boundary 
of r other than the Lipschitz condition, when T is an open surface. In any case, we believe 
that the restriction i/ > ^ is an artifact of our proof, and anticipating future weakenings 
of this restriction, the rest of this section will be presented so that it depends only on the 
assumption (35). 

Proof of Lemma 3.4- Let w € Sp satisfy 



u — w 



II 2 



+ osc(u;, P, rf < 2 inf (|||n - vf + osc(t;, P, Vf) . (40) 



veSp 
Then from the definition of oscillation, we have 

0SC{UP, P,Tf < 2 Y, h^r'lf - Mr,r + 2 Yl f^l'lM^ " ^p)lr,r 

T&P T&P (41) 

< 2osc{w, P,rf + 2Ca\Iw - upf , 

where in the last step we have used inverse inequality (35). This gives 

|||u - upf + osc(iip, P, rf < |||n - upf + 2Ca\Iw - upf + 2 osc{w, P, Vf, (42) 

and upon using the Galerkin orthogonality |||n — up|p + \w — up|p = |||u — wlp, we obtain 
(36) with, say Cq = 4(1 + Ca). 

Now we turn to b). With Aq G (0,1) being the contraction factor for h^ when r is 
refined once, we have 

E h"r\rp\lr< E >^th';\rp\l^ = \l^o.c{P'Xf. (43) 

T&P'\P T&P\P' 

Using this in 

osc(P', rf < (1 + <5) E K'\rp\lr + (1 + ^) E h?\A{up, - up)\l,, (44) 

tGP' t£P' 

and by using the inverse inequality (35) on the last term, we establish (37). 
For c), from (38) we infer 



(1 - 2/x)(|||ii - upf + osc(P, rf) < |||up - up>f + osc(P, F)2 - 2 osc(P', F)^ 

< (1 + 2Ca)\Iup - up>f + osc(P,F*)2 

< o \\rp\\r* +osc(P,F ) , 



(45) 
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where we have used the Galerkin orthogonahty, the estimate 

osc(p, r \ r*)2 < 2 osc(p', r \ r*)^ + 2Ca\Iup - up>f, (46) 

the local discrete upper bound in (33), and the norm equivalence (8). The proof is com- 
pleted upon noting that 

\\rpf <j\lu-upf, (47) 

where 7 = /3^/a > 1, which is a combination of the global lower bound in (13), and the 
norm equivalence (8). □ 

Once we have the preceding results, and given the techniques developed in Stevenson 
(2007) and Gascon et al. (2008), it is a fairly straightforward matter to obtain geometric 
error reduction and quasi-optimality for an adaptive method such as the ones considered 
here. Nevertheless, we include detailed proofs for convenience of the reader. Our first stop 
is the contraction property of the adaptive method. 

Proposition 3.6. Let the assumption (35) of Lemma 3.4 hold. Let P, P' G adm(Po) be 
adm,issihle partitions with P < P' , and let T* = int UrePVP' ^- Suppose, for some 9 £ (0, 1] 
that 

\\rp\\l, + osc(P, r*f > 6 {\\rpf + osc(P, r)^) . (48) 

Then there exist constants 7 > and fi £ (0, 1) such that 



u ■ 



np/|p + 7 0sc(P',r)2 <fi{\lu-upf + josc{P,Tf) . (49) 



Proof. The property (48), together with the global upper bound in (13), the local discrete 
lower bound in (33), and the norm equivalence (8), gives 

Q,2 \ 2/32 

'"n-upf + osc(P,r)M < -—\lup-up>f + 2{2Cj + lfosc{P,T*f. (50) 



Then for any 7 > 0, we combine Lemma 3.4 b) and (50) to infer 

III Ill2 I /J-)' T^\2 III |||2 III |||2 I /J-)/ T^\2 

\lu — upi\l +'yosc[P,L) =\lu — up\l — fup — up'f +'yosc[P,i) 

< ju - upf - (1 - -rCs)\lup - up,f + 7(1 + 5)osc(P, rf - 7A(1 + 6)osc{P, r*f 
( ea^{l-^Cs) \ 2^( t-,^,, ea{l-^Cs) \ .p^.2 

- V W^ — J ~ V W^ — J 



a(2Cj + 1)2(1 -7Q) w, ^..\ ,pp..2 
"2 7A(1 + 6)\ osc(P, r ) 

=: /xilfu - upf + ^l2l osc(P, Vf + /i3 osc(P, T*f. (51) 



3 Operators of order zero 15 



We choose 7, depending on 5 > 0, so that /X3 = 0, i.e., so that 

l^C.^^g^. (52) 

7 a{2Cj + ly 

With this choice and for 5 > sufficiently small, we have 

^' = '- 2i2Cj + ir <^' ^''^ 

and 

establishing the proof. □ 

To be explicit, the importance of the preceding proposition is the following. Suppose 
that we have an admissible partition P and the discrete solution up G Sp. Then the local 
residual norm HrpH,- and the local oscillation osc(P, r) are computable for t £ P, and 
by selecting a set i? C P of triangles such that T* = int Urs-R ^ satisfies (48) for some 
6 E (0, 1), and then by finding an admissible refinement P' of P such that P\P' D R, i.e., 
that each triangle in R is refined at least once, we can guarantee e(P') < fJ-e{P) for some 
/i < 1, where e{Q) = |||n — uqIP + josc{Q,T)'^. By repeatedly applying this procedure, 
we can ensure convergence |||m — up||| — )• as P runs over the partitions generated by the 
algorithm. This however, does not say anything about the growth of i^P, a fundamental 
question that will be addressed below. 

For u G L^(r) the solution of Au = f with / G H'^{T, Pq), and for P £ [Pq], we define 

dist^('u, Sp) = min (l-u - vf + osc(w, P, Vf) ^ , (55) 

v^Sp 

where the implicit dependence of oscillation on r has been made explicit. Note that the 
minimum exists since Sp is finite dimensional. Furthermore, for e > we define 

cardr.(ti,e) = min{#P - #Po : P G adm(Po), distr(u,5p) < e}. (56) 

Hence cardr.(u,e) is in certain sense the cardinality of a smallest admissible partition P 
that is able to support a function that is within an e distance from u. Note that cardr(n, e) 
is finite for any e > 0, since from the discussion in the preceding paragraph, there is a 
sequence of (finite) partitions {Pk} C adm(Po) with distr(u, Sp,.) — )■ 0. 

Another notion we need is that of overlay of partitions: We assume that there is an 
operation © : adm(Po) x adm(Po) — ^ adm(Po) satisfying 

P(BQhP, PeQhQ, and #(P © Q) < #P + #Q - #Po, (57) 
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for P,Q€z adm(Po)- In the conforming world, P (B Q is taken to be the smallest and 
common conforming refinement of P and Q, for which (57) is justified in Gascon et al. 
(2008). 

The following proposition shows an optimal way to select the triangles to be refined. 
The procedure is known as Dorfler's marking strategy in literature. 

Proposition 3.7. Let the assumption (35) of Lemma 3.4 hold. Let P G adm(Po); o,nd 

3 

let 9 G (0,^*) with 6* = „3^° ^ . . Suppose that R ^ P is a subset whose cardinality is 
minimal up to a constant factor k > 1, among all R C P satisfying 

\\rp\\hiR) + osc(P, r*{R)f > e {\\rpf + osc(P, Tf) . (58) 

where T*{R) = intjj^g^r. Then we have 

#R < Kcardr{u,e), (59) 

where e is defined by 

e' = ^^{\lu-upf + osc{P,rf). (60) 

Proof. Let us introduce the abbreviation 

e(Q) = |||'u-'UQ|||2 + osc(Q,r)2, Q e [Po], (61) 

and let P^ £ adm(Po) be such that 

#p^ - #Pq < cardr.(ii,e), and e{Pe) < e^ . (62) 

Then for P = P (B Pe, from Lemma 3.4 a) we have 

e{P)<CGe{Pe)<CGe^ = fie{P), (63) 

where fi = -i^^ G (0, ^), and so an application of Lemma 3.4 c) gives 

lkp|lr*(P\P) + os<^' r*(P \ P))2 > e (||rp||2 + osc(P, r)2) . (64) 

Recalling that R minimizes f^R up to the constant factor k among all subsets P \ P 
satisfying the preceding inequality, we infer that #i? < k#{P\P), and taking into account 
that #{P\P) <#P-#P and the estimate in (57), we get 

#R < k(#P - #P) < K{#Pe - #Po) < K card,(u, e), (65) 

completing the proof. □ 
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We have almost all the ingredients to give a bound on the growth of #P, and to discuss 
whether this growth rate is optimal in one or another sense. Let us start by making the 
term "adaptive BEM" precise. 

Algorithm 1: Adaptive BEM {t = 0) 
parameters: conforming partition Pq, and 9 £ [0, 1] 
output : Pk S adm(Po) and Uk G Sp^ for all k £ Nq 

1 for A; = 0, 1, . . . do 

Compute Uk G Sp^, as the Galerkin approximation of u from Sp^ ; 

Identify a minimal (up to a constant factor) set Rk C Pk of triangles satisfying 



4 



5 endfor 



\\rk\\h + osc{Pk, T*f > 9 {Wnf + osc(Pfc, r)^) , (66) 

where r^ = f — Auk and T* = int Urei? ^i 

Set Pk+i as an admissible refinement of Pk with Pk \ Pk+i 5 Rk] 



In order to make Line 4 in this algorithm unambiguous, we assume that a procedure 
to compute Pk+i € adm(Po) from Pk and Rk is specified, and that this procedure satisfies 

k-i 
#Pk - #Po <CcJ2 #^-' ^ ^ 1' (67) 

for any given sequence {Rm} subject to the obvious constraint Rm C Pm, where Cc > 
is a constant. This assumption is justified for newest vertex bisection algorithm in Binev 
et al. (2004); Stevenson (2008), and moreover, can be demonstrated for the red refinement 
procedure by adapting their techniques. 

Now we turn to the issue of convergence rate. Since Proposition 3.7 gives a bound on 
#i?fc, we simply need to use (67) to get 

k-l 

#Pk - #Po < kCc Y^ card,(n, Ce{Pm)), (68) 

m=0 

where C is the constant from (60), and e{Pm) is as in (61). On the other hand. Proposition 
3.6 guarantees a geometric decrease of e{Pm), i.e, we have 

e{Pk)<Cfi''~"'e{Pm), hence card^(ii, Ce(P„)) < card^(u, C^'"-^e(Pfc)), (69) 

for some constants C > and < /i < 1. Note that C denotes different constants in its 
different appearances. Therefore, if our particular u G A~^ {H^ {T , Pq)) satisfies 

card,.(M, Xe) < CX'^^^c&rdr{u, e), A > 1, e > 0, (70) 
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for some constant s > 0, then we would get what can be cahed instance optimality 

#Pk - #Po < Ccard,(n, e{Pk)). (71) 

However, it is not clear how to usefully characterize the set of such u, so we settle for 
something more modest. Instead of (70), if we have 

cardr(n,e) < C„e"^/^ e > 0, (72) 

for some constant s > 0, then we get what can be called class optimality 

#Pk - #Po < CCueiPkY^/'. (73) 

This motivates us to define the approximation class Ar^s C A~^ {W^ {V , Pq)) with s > 
to be the set of u for which 

|w|^^^ = sup (cardr(u, e)^e) < oo. (74) 

£>0 

Thus Ar,s is characterized by 

cardr.(n,e) < e"^/''|uU,,^^, (75) 

or equivalently, by 

min{dist^(u, 5p) : P e adm(Po), #P - #^0 < A^} < N~'\u\a,^^. (76) 

We have proved the following. 

Theorem 3.8. Let the assumption (35) of Lemma 3.4 hold, and in Algorithm 1, suppose 
that £ (0, 6*) with 9* = gww^2Cl- ^^^ ^ ^ H'^iT, Pq) and u G A,s for some s > 0. Then 
we have 

|||n-nfc|||2 + osc(nfc,Pfe,^)2<C|n|3^^J#Pfe-#Po)"'^ (77) 

where C > is a constant. 

Remark 3.9 (Approximation classes). The preceding theorem establishes optimal conver- 
gence rates for our adaptive BEM with respect to the class Ar,s- However, this approxima- 
tion class is based on a modified notion of error (55) that has the oscillation term in it, while 
what is really relevant is the class As that is based on dist(n, Sp) = m(y,^Sp fu — vf. Note 
that if r = 0, then the oscillation is dominated by the error, so As = Aq^s with equivalent 
(quasi-) norms. As far as the author is aware, the idea of a modified approximation class 
for adaptive methods first appeared in Gascon et al. (2008). This is an important idea as 
it allows for breaking up the optimality problem into two natural parts: first to establish 
optimality in some approximation class that is naturally tied to the particular algorithm, 
and then to show that those classes are sufficiently large so as to contain the corresponding 
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standard approximation classes. The approximation class Ag proposed by Gascon et al. 
(2008) can be adapted to our setting as the class of triplets {u,f,A) satisfying u G Ar,s, 
f = Au £ -ff''(r,Po)) and A varying in some predetermined class of boundary integral 
operators. There seem to be some convenience in using classes of singletons n, rather than 
triplets (u, /, A). First of all, / is redundant since / can be deduced from / = Au. More 
importantly, we want the approximation classes to be linear spaces as opposed to nonlin- 
ear sets of triplets {u,f,A) or pairs {u,A). Put differently, we view the nonlinear space 
As as a collection of linear spaces parameterized by A. Then we propose the problem of 
characterizing Ar,s on a large relevant class of A, which is a more tractable problem than 
giving a complete characterization of Ag. 

4 Positive order operators 

In this section, we consider the case t G [0, |). We assume that F is open, and remark 
that the closed case can be treated with similar methods. Recall that the equation we are 
dealing with is Au = f with linear homeomorphism A : H^{T) — )■ H~^{T). For simplicity we 
consider right hand sides satisfying / G L^(F). We employ the continuous piecewise affine 
functions with homogeneous boundary conditions Sp := Sp. The admissible partitions 
will be used synonymous with conforming triangulations, with the newest vertex bisection 
algorithm for refinements. So we have Sp C if^*(F), and we make an additional assumption 
that A : H'^^{T) — ;> i^(F) is bounded, in order to keep the useful property f — Av G i^(F) 
for any finite element function v G Sp. Note that this assumption is satisfied for our main 
example ~ the hypersingular integral operator. Note also that if t > | then we need u >2, 
i.e., we need the space fi in which F lies to be at least a C^'^ manifold. 

Starting from this section, we shall often dispense with giving explicit names to con- 
stants, and use the Vinogradov-style notation X <Y, which means X < C -Y with some 
constant C that is allowed to depend only on the operator A and on (geometry of) the 
set of admissible partitions adm(Po)- The following is an extension of Lemma 3.1 to f > 
and piecewise linear finite element spaces. In contrast to Lemma 3.1, note that co* includes 
a buffer layer of triangles around the refined triangles, so that uj* = a;(F*) with F* as 
in Lemma 3.1. By using the Scott-Zhang quasi-interpolation operator in the proof, it is 
possible to do without the buffer layer if t > ^ , which is however not terribly exciting since 
the relevant case to us is t = g- 

Lemma 4.1. Let P,P' G adm(i-o) "^ith P ^ P' , «^^ ^ei uj* = Utgpvp' '^(''")- ^^en we 
have 

lup - up>l < \\hprp\\^* < ||/ip//ip/||*^|||np - up>\l + ||/iprp - /ip*w||^*, (78) 

for any v G Sp' with suppu C n;*, where \\ ■ ||oo denotes the L°°-norm. Moreover, we have 
the following global bounds 



u ■ 



up\l < \\hprp\\ < \\hp/hp,\\l^\lu - upj + ||/iprp - hp^v\\, (79) 
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for any v G Spi . 

Remark 4.2. The global upper bound (i.e., the first inequality) of (79) has been established 
in Carstensen, Maischak, Praetorius, and Stephan (2004), and a similar bound involving 
local L^-norms on stars appears in Nochetto, von Petersdorff, and Zhang (2010). 

Proof of Lemma 4-1- For v £ -ff*(r) and vp = Qpv G Sp being the quasi-interpolant of v 
as in (22), we have 



{rp, v) = {rp, V -vp) = {hprp, hp\v - vp)) 

< \\hprp\\\\hp\v - vp)\\ < \\hprp\\\\v\\jjtf^^y 

where in the last step we have used 



(80) 



J2 K''\\v - vp\\l < Y, lbllH*(-M) ^ ll^lli(r)' (81) 

reP reP 

by local finiteness. This gives the first inequality in (79). 
With V = upi — Up and vp = Qpv £ Sp, we infer 

{Av, v) = {rp, V -Vp) = {rp, v - vp)^» 

< \\hprp\\^*\\hp\v - vp)\\:^* < ||/i*prp||c^. ||w||j^t(r), 



(82) 



where we have used the fact that vp = v outside uj* . This gives the first inequality in (78). 
Let either e = up^ — up oi e = u — up. Then in both cases, for v G Sp' and w D supp u, 
we have 



{hp^v, hp^v)^ = {hp^^v, v) = {hp^^v - rp, v) + {Ae, v) 
= {hp^v - hprp, hp^v) + {Ae, v) 

< ||/ip*f - hprp\\^\\hp'^v\\ui + \\Ae\\H-t(^r)\MHt{r) 

< ||/ip*u - hprp\\i^\\hp^v\\^ + ||^e||/^-t(r)||/ip*'v||a;, 



(83) 



so that 



II V^IU ^ II V^ - h^prpWuj + II V/V|lLDl|e||^t(r)- (84) 
From this, we infer 

\\hprp\\^ < \\hprp - h~p^v\\u; + \\hp'^v\\^ < \\hprp - hp^v\\^ + ||/ip//ip'||Ll|e||j^t(r), (85) 

proving the both second inequalities in (78) and (79). □ 
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Remark 4.3 (Operators of order zero). The preceding result has important impUcations 
also in the case t = 0. Assuming that / G -f^^(r, Pq), we can define the oscillation 



osc(P,a;)= I >^ hl\rp\l\ , (86) 




for a; C r, so that (78) and (79) imply 

jup-up'W < \\rp\\^* < jup-up'W +osc(P,tj*). (87) 

and 

III^-^pIII < Ikpll < lll^i-wp|||+osc(p,r). (88) 

Then the techniques from Section 3 can be adapted to obtain a quasi-optimal convergence 
rate in this setting as well. The only new features are the oscillation term in the global 
upper bound, and the buffer layer of triangles in uj* as compared to T* . 

Remark 4.4 (Saturation assumption) . For a large number of non-residual a posteriori error 
estimators for the hypersingular integral equation on curves, Erath, Funken, Goldenits, and 
Praetorius (2009b) proved that the estimators are equivalent to the global error, with the 
upper bound depending on the saturation assumption: |||ti — upf < |||np — upf, where up 
is the Galerkin approximation from some enriched space Sp D Sp, which is typically the 
piecewise linears on the uniform refinement of P. Combining the discrete lower bound with 
the global upper bound from Lemma 4.1, we have 

|||u — Up III < III Up — Up/ III + ||/iprp — hp^v\\, (89) 

for any v G Sp', where P' is the uniform refinement of P. This confirms the saturation 
assumption up to an oscillation term. Then adapting the techniques from Ferraz-Leite, 
Ortner, and Praetorius (2010), one can expect to get a geometric error reduction for adap- 
tive algorithms based on non-residual type estimators; however, it seems that any type of 
quasi-optimality result would require new ideas. 

Let us get back to the case t > and the residual based error indicators || /i^rpH,- from 
Lemma 4.1. We define the oscillation for v £ Sp by 

osc{v,P) = ||/i*p(/ - Av) - hp^wl (90) 

where, with P the uniform refinement oi P, w = Qph?p{f — Av) E Sp is the Clement 
interpolator of /ip(/ — Av), given by 

Qp9= Yl 9z{z)(t>z, (91) 

z&Np\dV 
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where Np is the set of all nodes in P, (j)^ G Sp is the standard nodal basis function at z, 
and Qz G IPi is the L^-orthogonal projection of g onto the affine functions on the star C^{z) 
around z with respect to P. Then for 7 C F and 1^(7) = lj{d)(z) : z G Np Pi 7}, we have 

Wh'prp - hp'Qphj^rpf^ < CgWh'prpWl^^y (92) 

for some constant Cq > 0, i.e., the estimator dominates the oscillation on the Galerkin 
solutions. In particular, we have the equivalence 

ai\\h'prp\\^ < |||n-upf + osc(up,P)2 < \\h'prpf, (93) 

where ai > is a contstant. Hence the situation is very similar to what happens in the 
finite element case, and in particular, the local quantities || /ipJ^pUr as error indicators will 
give rise to an adaptive algorithm that converges quasi-optimally in a certain sense. 

Lemma 4.5. Assume that 

Wh^pAvf < CAhf, V e Sp, (94) 

for P G adm(i-o)) with the constant Ca = CAiA,as,(Jg). Then the followings hold. 

a) There is a constant Co > such that 

ju-upf + osc{up,Pf <Cg ini {ju - vf + osc{v , Pf) , (95) 

v&Sp 
for any P £ adm(Po)- 

b) There exists a constant A > such that 

\\h'p,rp,f < (1 + 6)\\h'prpf - A(l + S)\\h'prp\\l, + Csjup - up4'^, (96) 

for any P,P' G adm(Po) with P :< P' , and for any 5 > Q, with Cg depending on S, 
where T* = int UrePXP' ^• 

c) Let P,P' G adm(Po) be such that P :< P' , and let T* = int |J^gp\ p, r. If, for some 
/i G (0, 2) it holds that 

\u — upi\ + osc{up', P') < /i (III If — Up III +osc(up,P) ), (97) 

then with 6 = ^ j^ 11I2C C ) ^^^ ^* ~ UreP\P' '^(''")' "^^ have 

\\h'prp\\l,>e\\h'prpf, (98) 

where f3i > is a constant implicit in the local discrete upper bound in (78), such that 
\lup-up>f<Pi\\h'prp\\l,. 
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Remark 4.6. The estimate (94) is proved in Section 6 for a general class of singular integral 
operators of positive order, under the hypothesis that A : H''^^(il) — t- H^*^'^{Q) is bounded 
for some a > max{^,t}, where A is an extension of A to fl. Recalling that Q is a C'^~^'^- 
manifold, and assuming that t < ^, the condition on a translates to z^ > § + i, see 
Remark 6.3 for details. Therefore for n = 1 and n = 2 we need a C^'^ curve or a surface, 
respectively, assuming that ^ < t < 1. Note that the boundary of F is allowed to be 
Lipschitz. Anticipating future developments on (94), the rest of this section is presented 
so as to depend only on the assumption (94). 

Proof of Lemma 4-5. For v G Sp we have 

osc{up,P) = \\hprp - hp^QphprpW 

< osc{v,P) + \\h^pA{v - up) - h^^Qph^^A{v - up)\\ (99) 

< osc{v,P) + C]^'^\\h^pA{v - up)\\ < osc{v,P) + &J^C]/^\lv - up|||. 

With this estimate at hand, part a) of the lemma can be proven in exactly the same way 
as Lemma 3.4 a). 

Part b) follows from 

\\hp,rp>f < {l + 6)\\h^p,rpf + CA{l + S-^)\lup-up'f, (100) 

and the fact that hp/ < Xohp on F* for some constant Aq < 1. 

For c), the global lower bound, property (97), and the Galerkin orthogonality give 

ai(l - 2;u) II /i^rp II 2 < (1 - 2/x) {ju - upf + osc(np, Pf) 

< l\u - upf + osc{up,Pf - 2|||m - up'f - 2osc{up',P'f (101) 

< |||up — Mp'l +osc(mp,P) — 2osc(mp', P') . 

For T £ P f] P' with all its neighbors also in POP', where P' is the uniform refinement 
of P' , we have {Qpv)\r = {Qp/v)\t, and also hp = hpi there. Hence, with cb = (J{r G P : 
3<T G P \ P', T n CT / 0}, we have 

||/iprp - /ip*(5p/ip?'p||r\<i < \\hpirp - hp^Qp,hf,rp\\ 

< \\hp,rp' — hjJ;Qp,hp,rp'\\ + \\hp,A{up — up/) — hjJ;Qphp,A{up — up')\\ 

< osc(tip/,P') + Cq''^||/i^,A(up - up/)||. (102) 

Combining this with (92), we infer 

osc(np,P)2 < 2osc(up/,P')^ + CQ||/i*prp||2, +2CQC7A|||up-np/f , (103) 

which, then on account of (101) and the local discrete upper bound, proves the claim. □ 
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Now we prove an analogue of Proposition 3.6 on error reduction. 

Proposition 4.7. Let the assumption (94) of Lemma 4-5 hold. Let P,P' G adm(Po) be 
with P :< P' , and let V* = int |JreP\p' ^- Suppose, for some 7? G (0, 1) that 

||/ikp||r*>^||/ikp|l- (104) 

Then there exist constants 7 > and fi £ (0, 1) such that 

|||'u-Mp/|P + 7||/i*p,rp/|p < /u(|||'u-'up|p + 7||/i*prpf) . (105) 

Proof. From the Galerkin orthogonality and Lemma 4.5 b) we have 

III O 11+ 1 1 'J O III 'J 1 1 ■/■ 1 1 

|m — np/| + 7||/ip'?^p'|| =iF~'Wpi ~iFP~'"P'i + tII^p'^'p'II 

< |||u - upf + 7(1 + 6)\\h'prpf - 7A(1 + d)\\h'prp\\l, + {-fCs - l)|||np - up.f 

< |||n - upf + 7(1 + 6)\\h'prpf - 7A(1 + (5)||/i^rp||2., (106) 

for < 7 < 1/Cs. The idea, introduced in Gascon, Kreuzer, Nochetto, and Siebert (2008), 
is to use the global upper bound and the Dorfler property (104) to bound fractions of the 
first two terms by the third term. For any a, 6 > be have 

|||m — lip'l + 7||/ip'?^p'|| 

< (1 - -/a)\lu - upf + 7(C7a + 1 + 5)\\h'prpf - -fX{l + 6)\\h'prp\\l, 
< (1 - 70)111^ - upf + j{Ca + 1 + 6- b)\\h'prpf 

+ 7(^-A(l + 5))||/ikp|lr., (107) 

with the constant C coming from the global upper bound. We choose b = A??^ > so that 
the third term is negative no matter how small (5 > is, and then choose a > and 6 > 
so small that Ca + 6 < b. □ 

For u G H^(r) the solution of Au = f with / G L^(r), and for P G adm(Po)) we define 

disto('u, Sp) = inf (Ik - vf + osc(w, Pf) ^ . (108) 

veSp ' 

Furthermore, for e > we define 

cardo(u,e) = min{#P - #Po : P G adm(Po), disto(n, S'p) < e}. (109) 

The following is an analogue of Proposition 3.7. 
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Proposition 4.8. Let the assumption (94) of Lemma 4-5 hold. Let P G adm(Po); o.iT'd let 
£ (0, 9*) with 6* = ^ g (il_2C c ) • Suppose that R^ P is a subset whose cardinality is 
minimal up to a constant factor, among all R C P satisfying 

\Wprp\\l,^^^>e\\h'prpf, (110) 

where T*{R) = intlj^g^r. Then we have 

#fl<cardo(n,e), (111) 

where e is defined by 

e' = ^^ (in - upf + osc(P, r)^) . (112) 

Proof. Although this proof is almost a copy of the proof of Proposition 3.7, there is a 
subtlety caused by the buffer layer of triangles in uj* as compared to T* . Let 

e{Q) = \lu-UQf + osc{Q,rf, Q G adm(Po), (113) 

and let Ps G adm(Po) be such that 

#Pe - #Po < card,(u,e), and e{Pe) < e\ (114) 

Then for P = P (B Pe, from Lemma 4.5 a) we have 

e{P)<CGe{Pe)<CGe^ = fie{P), (115) 

where fi = -^^ G (0, 5), and so an application of Lemma 4.5 c) gives 

\\rp\\l,>e\\rpf, (116) 

where u* = (J^gpvpa;(T) = T*{R), with R being the set of triangles in P who has a 

neighbour in P \ P. Since R minimizes ^R up to a constant factor, among all subsets 
R <Z P satisfying the preceding inequality, we infer 

#R<#R< *{P \P)<#P-#P< #Ps - #Po < card,(n, e), (117) 

where we have used the local finiteness and the estimate in (57), completing the proof. □ 

For completeness, in the rest of this section we give an explicit pseudocode for the 
adaptive BEM, then define the relevant approximation classes, and finally record a theorem 
on quasi-optimality. 



2 
3 
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Algorithm 2: Adaptive BEM (t > 0) 
parameters: conforming partition Pq, and 9 G [0, 1] 
output : Pk G adm(Po) and Uk € Sp,. for all fc G No 

1 for A; = 0, 1, . . . do 

Compute Uk G Sp,, as the Galerkin approximation of u from 5p^ ; 

Identify a minimal (up to a constant factor) set R^ C P^ of triangles satisfying 

Wh'p^TkWr* > 0\\h'pjk\\, (118) 

where r^ = f — Auk and F* = int UtS-R ^' 

Set Pfc+i as an admissible refinement of Pk with Pk \ Pk+i 5 Rk] 
5 endfor 

As in the previous section, we assume that a procedure to compute Pk+i G adm(Po) 
from Pk and Rk is specified, satisfying the bound (67). This is justified for newest vertex 
bisection algorithm in Binev, Dahmen, and DeVore (2004); Stevenson (2008). We define 
the approximation class Aq^s C A~^{L'^(r)) with s > to be the set of u for which 

\u\_Ag s = sup (cardo(n, eYe) < oo. (119) 

e>0 

The following result is immediate. 

Theorem 4.9. Let the assumption (94) of Lemma 4-5 hold, and in Algorithm 2, suppose 
that 6 G (0,6**) with 6* = ^ _^_„ (1I2C C ) ■ ^^^ ^ ^ -^^(r) '^'"■d u G Aq^s for some s > 0. 



Then we have 

IIL. „, ll|2 I „„„/„. D ^2 ^ /^L., 

where C > is a constant. 



u - Ukf + osc(ufc, Pkf < C\u\l(#Pk - #i^o)~'^ (120) 
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Finally, we turn to the case t < 0. The domain F can be either a closed manifold or a 
connected polygonal subset of a closed manifold. Recall that we are dealing with the linear 
homeomorphism A : -ff*(F) — t- H~*{r). We will use the piecewise constant finite element 
spaces Sp := Sp, and take conforming triangulations as the class of admissible partitions. 
Faermann (2000, 2002) established the equivalence 

Ilk - upf < 5^ |rp|2_,,,(,) < |||u - upf, (121) 



zeN 



p 



where Np is the set of all vertices in the triangulation P, and proposed to use the local 
quantities |rp|_j,^(^) as error indicators for adaptive refinements. An alternative proof, us- 
ing interpolation spaces appeared in Carstensen, Maischak, and Stephan (2001). However, 
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the discrete local counterparts to this equivalence have been open. On the other hand, 
the weighted residual type error indicators, h].~^^\rp\i^r for t £ P, have been around for a 
while, with a guaranteed global upper bound 

\lu-upf<^hl(^^^)\rp\l^, (122) 

cf. Carstensen et al. (2001). These indicators are computationally more attractive, but 
for locally refined meshes no global lower bound or discrete local estimates seem to be 
known. There is a recent claim though on two of those estimates, in Feischl, Karkulik, 
Melenk, and Praetorius (2011). The following lemma establishes the missing bounds for 
the afore-mentioned error indicators. 

Lemma 5.1. Let P,P' G adm(Po) ^e with P < P' , and let V* = intU^gp\p, r. Further- 
more, with uj* = uj{T*), let (p £ Sp, be the function such that (j) = 1 on T* and = 
outside oj* . Then for any v £ Sp,, it holds that 

|||Mp-up/f < Y^ \rp\'it^u;(,) + \\hp,irp-v)\\l* + \\(j)rp-vfj^^t(^r)- (123) 

zeNpHu)* 



In addition to the above, for some r G [— t, |) n (0, ^ — 2t) n [0,z^], let f G H^'{T) and let 

i. Then 
uniform refinement of P, we have 



A : if''+^*(r) — )• H^{r) be bounded. Then for any v G Sp and w G Spi, where P is the 



zeNpHu)* zGNpHu)* 

+ \\h'pirp-w)\\l + \\hUrp-v)\\l,+ Y. hf'^'''^\^p-v\l^^^^, (124) 

zeNpnoj* 

where h^ = mini^.^g;^! hr- Moreover, for any v G S^~, we have the global bounds 

zeNp zeNp 

+ Wh'pirp -v)f+Y hl^'''''^\rp - v\Uzy (125) 
zeNp 

Remark 5.2. In a previous version of this article, the local discrete upper bound (123) in- 
volved the factor \\hpi /hp\\'^ on the left hand side (which allowed to remove the oscillation 
terms) . It does not cause trouble as long as one is concerned with refinements P' of P that 
are not too deep, such as the ones that occur in practice. However, it is not effective in the 
current approach to optimality proofs as one has to consider arbitrary refinements of P, cf. 
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Lemma 3.4 c), Lemma 4.5 c), and Lemma 5.6 c) below. This observation was communi- 
cated to the author by Dirk Praetorius. In the current version, the afore-mentioned factor 
is removed with the help of the oscillation terms, and the rest of this section is rewritten 
to reflect this change. 

Proof of Lemma 5.1. Set s = — t > 0. Let e = upi — up, and let ep S Sp be the L^- 
orthogonal projection of e onto Sp. Then for any v G Sp,, we have 

{Ae, e) = {rp, e) = {rp, e-ep) = {rp, e)r* = {rp - v, e)r* + {v, e)r* 
< \\hpf{rp - v)||r*||/ip'e||r* + Ibllfl-^llellj^-s, 

where we have used the fact that e — ep = outside T* . Upon using the inverse inequality 
ll^p'^ll ^ IkllH-s; this gives 

ll|e||| ^ \\hpf{rp - v)\\r* + \\v\\hs < \\hpf{rp - v)\\t* + Urp - v\\hs + UrpWns. (127) 

To localize the last term, we follow an approach from Carstensen et al. (2001). With 
(pz G Sp the standard nodal basis function at z £ Np, we have 

\\(f>rp\\H- < Y^ Uzrp\\Hs(c(z)), (128) 

zGNpnuj* 

and then using the fact that the linear operator T : H^(uj{z)) — )• H^{uj{z)) defined by 
Tf = {f - (/, 1)^(^))(A2 satisfies \\Tf\\Hs < \f\H^uj{z)), we get the upper bound (123). 

For the lower bounds, we first prove a couple of inequalities involving the auxiliary 
error indicator ||/ip*rp||,-. Let w S Spi, and let wp S Sp be the L^-orthogonal projection 



of hp ^w onto Sp. Then we have 



{hp^w, hp^w) = {w — rp, hp ^w) + {A{up' — up), hp ^w — wp) 

< {w-rp,hp'^''w) + \\A{up/ - up)\\Hs(r*)\\hp'^''w - wp\\fj-s(^p,-) (129) 

< Whp^^w - rp)||suppt«||/ip''u'|| + \\up/ - up\\H-=(r)\\hp''w\\r*, 
where in the last step we used (20), and so 

Whp'^rpWsuppw < jup - up'W + \\hp%rp - tt;)||supp«,- (130) 

Let V £ L'^, and let vp E Sp be the L^-orthogonal projection of v onto Sp. Using (20), 
then we have 

{hr^rp, h^v) = {rp, v) = {A{u — up),v — vp) 

< \\A{u - up)\\hs{t)\\v - vp\\fj.s(Y) < |||n-'up|||||/i''v||, 

establishing 

\\h~p'rp\\<\\u-up\\. (132) 
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On the other hand, for any N C Np and v G S^, we have 

< 'S^ /"ll/j^^ilP _L?,2{r-s)| |2 \ 

~ Z^ l^ll'^P ^IL(^) +"2 FP ^lr,aj(2)j 

zGAf 



z^N z&N 



(133) 



zeAT 



< \\h-/rpf^ + \\hp%v - rp)||2 + ^ /i2(-^)|rp 



'^lr,a;(2)' 



z£N 



where 7 = int IJzgAr w(z). Then the second inequahties in (124) and (125) follow from (130) 
and (132), respectively. Finally, the first inequalities in (124) and (125) are a consequence 
of the fact that rp is L^-orthogonal to Sp. □ 

Let us record some useful bounds on the various oscillation terms that appeared in the 
preceding lemma. 

Lemma 5.3. Let P,P' £ adm(Po) be with P :< P' , and let 7 = intlj^ggr with some 
Q <^ P. Assume that f E H''{T) and that A : H''+^\T) -^ H''{T) is bounded for some 
r £ (-i, §) n (0, 5 - 2t) n [0,iy], so that rp £ H''(T). Then we have 

min \\h'p{rp - w)f^ < J^ Z^^^^+^VpI'.- (134) 

With V = Qprp £ S^~ the quasi-interpolant of rp, we also have 

zeNpH'j zeNpCrj 

where h^ = mini^.^g:^! h^. Finally, with T* and u* as in Lemma 5.1, there exists v £ Sp, 
such that V = outside uj* and that 

\\hUrp-v)\\h + Urp-v\\l^,^^^< Yl ^f'^'^^l^PlU^Y (136) 

zeNpHuj* 

Proof. The estimate (134) is a standard direct estimate, and (135) follows from 

E \\rp - v\\\^^.) < E hf'^'^\rp\U.iiz))< E hl^'^'^\rp\Uz), (137) 

zeNpn'j zeNpH'y zeNpCrj 

where we have used (21), (17), and the local finiteness of the mesh. 
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For (136), let v G Sp, be defined by 

v= Yl {'rp,l)uj'{z)4>z, (138) 

zeNp/Duj* 

where (j)z G Sp, is the standard nodal basis function at z G Np', and (•, ■)uj'(z) is the L'^-inner 
product on the star uj'{z) around z with respect to P' . Then we have 

\\hUrp-v)\\h< Yl \^pW^'(z)^ E \^P\-tMz), (139) 

zeNpiCiuj* zeNpDuj* 

and from rp _L^2 Sp we infer the bound (136) for its first term. For the second term, we 
have 

(l)rp-v= Yl {rp-{rp,l)u.'{z))(t>z, (140) 

zeNpiDuj* 

and hence 

zeNp,nuj* 

As in the proof of Lemma 5.1, now by using the boundedness of / i— )• (/ — (/, l)aj'{z)) 4'z in 
H~*{uj'{z)), and then by employing the orthogonality rp ±1^2 Sp again, we establish the 
proof. □ 

Remark 5.4 (Saturation assumption). For a large number of non-residual a posteriori er- 
ror estimators for Symm's integral equation on curves, Erath, Ferraz-Leite, Funken, and 
Praetorius (2009a) proved that the estimators are equivalent to the global error, with the 
upper bound depending on the saturation assumption: fu — upf < |||mp — np|||, where up 
is the Galerkin approximation from some enriched space Sp D Sp, which is typically the 
piecewise constants on a uniform refinement of P, or piecewise linears on P. By using 
Lemma 5.1, we can justify this assumption up to oscillation when Sp is the space of piece- 
wise constants on a uniform refinement of P. Namely, the local lower bound (124) with 
r = —t, and the global upper bound from (121), together give 

h-upf<\lup'-upf + \\h^p{rp-w)f + \\h^p{rp-v)f+ Y W^P - ^W-tMz)^ (142) 

zeNp 

for any w G Sp' and v G S^~. Taking P' equal to P (the uniform refinement of P) implies 
the saturation assumption up to the oscillation terms. 

Remark 5.5 (Weighted residual indicators). Putting r = 1 in (124) and (125), one gets 
bounds for the weighted residual indicators from (122). Taking into account the preceding 
lemma, the oscillation terms in the global lower bound (125) can be estimated as 

WhUrp -v)f + Y h'r^'-'^l^P - -II. < Y hl^'-'^\^p\lr, (143) 

tGP tGP 
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for V £ S^~ minimizing the left hand side. Hence the estimator dominates the oscillation. 
This observation suggests that arguments similar to the ones given in the preceding section 
might lead to a proof of quasi-optimality of an adaptive algorithm based on the weighted 
residual indicators. A proof of a certain optimality result in this setting has been announced 
in Feischl et al. (2011), with the details appearing in a forthcoming paper. 

Combining (123) and (136), for the Faermann indicators we get 

zeNpHuj* zGNpHu)* 

By (124), the entire right hand side is controlled by the last term alone. This leads to 
generalized weighted residual indicators, for which we already have the global bounds (125). 
Obviously, the case r = 1, treated in Feischl et al. (2011) is computationally more attractive, 
but the remaining cases r G (— t, 1) become important if for instance / G H^(T) \ H^(T). 
In the following, we fix some r G {—t, |)n(0, ^ — 2i)n[0, u] and assume that / G H^{T), 
and that A : H^~^'^^{r) — t- H^{T) is bounded. Then defining the error estimator 



\zeNpny J 

for 7 C r and v G Sp, with r](P, 7) = r]{up, P, 7), in the context of Lemma 5.1 we have 

|||np-up'|||'</3i'?(i^,r*)2, (146) 

where /3i > is a constant we will refer to later. Let there be a map z 1— >• Tp{z) : Np — t- P 
for any P G adm(Po)) satisfying z G Tp{z). For v G Sp, let 

osc{v,P)= l\\h'p(f-Av-w)f+ E h'^-^'^lf-Av-wll^^J , (147) 

V z&Np J 

where co{z) is the star around z G Np with respect to the partition P, and w = Qp{f — 
Av) G S^~ is the quasi-interpolator of / — Av, defined by 

Qp9= E 9z{z)cl),, (148) 

zGAfp 

(pz G Sp is the standard nodal basis function at z, and Qz G Pi is the L'^-orthogonal 
projection of g onto the afhne functions on Tp{z). This quasi-interpolator, as a variation 
on the Clement interpolator, is introduced in Oswald (1994). We remark that we will 
only need the idempotence Q^~ = Qp, so for example the Scott-Zhang operator could have 
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been employed instead. Let g' = g — Qpg with g G H^{T), and let y £ Np. Then by 
idempotence, we have 

Wllu^iy) = la' - Qpg'lrMy) ^ Yl lid' - g'ziz))<Pz\rM^), (149) 

zeNpnuj{z) 
which implies by the boundedness of g' i— t- (g' — g'^{z))(j)z in H'''{lxj{z)), that 

E l5 - Qp9\l^{y) ^ E l5 - Qpa^'rMzy (150) 

y&Np z&Np 

In particular, in the context of Lemma 5.1 we have 

III'" - upf < airiiP, Tf < |||n - upf + osc{up, Pf, (151) 

where ai > is a constant. On the other hand, by stability of Qp there exists a constant 
Cq > such that 

iKirp - Qprp)\\l^,) + E_^'^'^*^l"^ - Qp^P\rMz) < CQViPnf, (152) 

zeNpnuj{j) 

for any P G adm(Po) and 7 ^ T, with ^(7) = IJ{'^(^) ■ ^ ^ -^p ^ 7}- 

In the rest of this section, we follow the pattern of the preceding two sections, and the 
proofs will be brief as similar proofs are given in the previous sections. 

Lemma 5.6. Assume that 

^ hl'^^+MvlU.) < CAhf, V E Sp, (153) 

zeNp 

for P G adm(Po)) with the constant Ca = ^^^(A, 0"^, 0"^). Then the followings hold. 

a) There is a constant Cq > such that 

ju-upf + osc{up,Pf <Cg ini (ju - vf + osc{v , Pf) , (154) 

v&Sp 
for any P G adm(Po)- 

b) There exists a constant A > such that 

viP', Tf < (1 + 6)viP, Tf - A(l + 6UP, T*f + Csjup - up,f, (155) 

for any P,P' £ adm(Po) with P :< P' , and for any 5 > Q, with Cs depending on 5, 
where T* = int UreP\P' ^- 
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c) Let P^P' £ adm(Po) be such that P :< P' , and let T* = int IJ^^pv p, r. //, for some 
/i G (0, 2) ii holds that 

\u — upi\ + osc{up',P') < /x (|||m — up||| + osc{up,P) ), (156) 

then with 9 = ^ — , a^i ~n^ r^ ^ , we have 

v{P,r*f>ev{P,rf, (157) 

where ai > and /3i > are the constant from (151) and (146), respectively. 

Remark 5.7. The estimate (153) is proved in Section 6 for a general class of singular 
integral operators of negative order, under the hypothesis that A : H^~^'^{Vt) — )• H~^~^"{Q) 
is bounded for some a > ^, and that < r + t < cr, where A is an extension of ^ to $7 if 
r is open, and A = AiiT = Q.. Recalling that il is a C"^~^'^-manifold, the condition on a 
translates to z^ > f — i, see Remark 6.3 for details. Therefore for n = 1 and n = 2 we need 
a C^'^ curve or a surface, respectively, assuming that ^ < — t < 1- Note that if V is open, 
the boundary of T is allowed to be Lipschitz. 

Proof of Lemma 5.6. The proofs of a) and c) go along the same lines as the corresponding 
proofs from the previous section, cf. Lemma 4.5. In particular, the stability (152) and the 
locality of Qp are important. 

Claim b) is established if we show that 

E_^^^^^*^I^IU)^^ Y. _hf^'^\9\Uyy 5^^^(r)> (158) 

with some constant A < 1, where ijo'{z) is the star around z G Npi with respect to the 
partition P' . For T,a £ P' , let I{T,a) denote the interaction term between r and a in 
the Slobodeckij (double integral) norm of g. First of all, any diagonal term /(r, r) that 
appears in the left hand side also appears in the right hand side, and the corresponding 
factors satisfy h^ < Xohy with some constant Aq < 1. Henceforth we concentrate on the 
off-diagonal terms. Note that by symmetry the order of r and a is not important, and 
that the number of occurrences of the particular (unordered) pair (r, a) in the left hand 
side of (158) is equal to the number of z G Np/ n F* satisfying r, cr C ijj'{z). Suppose that 
the pair (r, a) appears in the left hand side exactly i times, for I £ [0, n — 1]. Thus r and 
a are contained in two triangles from P that share a A:-face for some k £ [i — l,n], where 
k = n means that the two triangles coincide. If this face is in F*, then the vertices of this 
face give at least I points y G Np n F* such that r, a C t<j(y), meaning that the same pair 
appears in the right hand side at least i times. On the other hand, if the shared face is not 
in F*, this would mean that r and a are triangles from P, and they interact through only 
the vertices on Np n dV* . We also see that the corresponding factors hz shrink since 

hz is defined by taking minimum as hz = minr^gp/.^g^i /i,-, proving the claim (158). □ 
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Proposition 5.8. Let the assumption (153) of Lemma 5.6 hold. Let P,P' G adm(Po) be 
adm,issible partitions with P :< P' , and let T* = int Utgpvp' ^- Suppose, for some t!) G (0, 1) 
that 

ri{P,T*)>^ri{P,T). (159) 

Then there exist constants 7 > and fi £ (0, 1) such that 

|||n - up'f +-fv{P', r)^ < fi (|||n - upf +jr]{P,rf) . (160) 

Gearing towards a convergence rate analysis, for u G H^(r) the solution of Au = f 
with / G H^(T), and for P G adm(Po)) we define 

distriu, Sp)= inf (|||u - uf + osc(v, P)^)^ . (161) 

fSSp 

Furthermore, for e > we define 

card,.(u,e) = min{#P - #Po : P G adm(Po), distr{u, Sp) < e}. (162) 

We have the following result on the Dorfler marking, whose proof is entirely analogous to 
the proof of Proposition 4.8. 

Proposition 5.9. Let the assumption (153) of Lemma 5.6 hold. Let P G adm(Po)7 o,nd 
let 9 G (0,9*) with 9* = ^,0 (i\_2c ) ■ Suppose that R <^ P is a subset whose cardinality is 
minimal up to a constant factor, among all R^ P satisfying 

r]{P,r*{R)f>9r]{P,rf. (163) 

where T*{R) = intlj^g^r. Then we have 

#R<cardriu,e), (164) 

where e is defined by 



e^ = ^^^(\\\u- 



2Cg9* 
Now let us specify our adaptive algorithm. 



'Up|P + osc(P,r)2). (165) 



2 
3 
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Algorithm 3: Adaptive BEM {t < 0) 
parameters: conforming partition Pq, and 9 G [0, 1] 
output : Pk G adm(Po) and Uk G Sp,. for all fc G No 

1 for A; = 0, 1, . . . do 

Compute Uk G Sp,, as the Galerkin approximation of u from 5p^ ; 

Identify a minimal (up to a constant factor) set R^ C Pk of triangles satisfying 

rj{Pk,r*f>erj{Pk,rf, (lee) 

where rk = f - Auk and T* = int UreiJ^ "^5 

Set -Pfc+i as an admissible refinement of Pk with Pk \ Pk+i 5 Rk] 
5 endfor 

As in the previous sections, we assume that a procedure to compute Pk+i G adm(Po) 
from Pk and Rk is specified, satisfying the bound (67). This is justified for newest vertex 
bisection algorithm in Binev, Dahmen, and DeVore (2004); Stevenson (2008), and moreover, 
can be demonstrated for the red refinement procedure by adapting their techniques. Finally, 
we introduce the approximation class Ar,s C A~^ {H^ (T)) with s > to be the set of u for 
which 

|w|yi^^ = sup (cardr(M, e)^e) < oo, (167) 

£>0 

and record that our adaptive BEM produces optimally converging approximations. 

Theorem 5.10. Let the assumption (153) of Lemm,a 5.6 hold, and in Algorithm, 3, suppose 
that e G (0,^*) with 9* = 1+^^(1+20^) - ^e^ / ^ ^''(r) and u G A,s for some s > 0. Then 
we have 

\lu-Ukf + osc{Pk,Tf<C\u\X^^{#Pk-#Po)-^', (168) 

where C > is a constant. 

6 Inverse-type inequalities 

In this section we shall justify the inverse- type inequality (25), which has been used in 
Lemmata 3.4, 4.5, and 5.6, and hence played a crucial role in our analysis. We allow a 
general class of singular integral operators, specified by the assumptions that follow. 

We keep the assumptions formulated in Section 2 still in force. We will be concerned 
only with closed manifolds (i.e., F = 0), since the case of open surfaces F C would follow 
by restriction. In addition, we assume that Q is embedded in some Euclidean space M^, so 
that the Euclidean distance function dist : Jl x Q — )• [0, 00) is well defined. Instead of the 
operator A that featured in the previous sections, we will consider in this section a more 
general bounded linear operator T : H'^[Q) — )• //~*(S7), hence removing the self-adjointness 
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and coerciveness assumptions. With A = {(x, x) : x S 0} the diagonal of i7 x 0, we assume 
that there is a kernel K G L\^^{U, x $7 \ A) associated to T, meaning that 

{Tu,v) = {K,u®v), (169) 

whenever u,v ^ C"^~^'^(J7) have disjoint supports. We assume that K is smooth on 

S = {r X r' : T, r' e Po} \ A C J7 X f^ \ A, (170) 

satisfying the estimate 

for all multi-indices q and /3 satisfying ?7- + 2i + |a| + |/3| > 0. Note that the partial 
derivatives are understood in local coordinates (or, as we discussed in §2, in terms of 
the reference triangles). The kernels satisfying this smoothness condition has been called 
standard kernels, e.g., in Dahmen, Harbrecht, and Schneider (2006). Then one can show 
that the kernels of a wide range of boundary integral operators are standard kernels, cf. 
Schneider (1998). 

Theorem 6.1. With T* denoting the adjoint ofT, let both T,T* : H^+''{Vl) -^ F"*+'^(il) 
he hounded for some a > ^. Moreover, let s >0 and < s + t < a. Then we have 

j; /i^(^+*)||T^||^.(^(,)) < Ml.^^y V e Sp, (172) 

zeNp 

for P G adm(Po); where Np is the set of all vertices in the triangulation P, oo{z) is the 
star around z with respect to P, and h^ = minr^.^g^j hr- 

Corollary 6.2. Under the setting of the preceding theorem, we have 

J2hr^'^'^\\T^\\h(r) < Mlm^ V e Sp, (173) 

TdP 



for P G adm(Po)- 

Remark 6.3. The boundedness of both T,T* : H^^"{Vl) — t- H^^'^^{^) has been proved for 
o" < 2 for general boundary integral operators on Lipschitz surfaces in Costabel (1988), 
and the endpoint case a = ^ is established for boundary integral operators associated to 
the Laplace operator on Lipschitz domains in Verchota (1984). Unfortunately, we see that 
the preceding results require more than o" = g- Since v > a + |t|, we necessarily have 
u > ^ + |t|. This allows Lipschitz curves for |t| < ^i and C^'^ curves and surfaces for 
\t\ < 1. Even though in general it rules out polyhedral surfaces, note that for the case of an 
open surface T, its boundary can be Lipschitz polygonal, as long as one can find a smooth 
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enough manifold O with F C il. As for the question of whether the boundedness holds 
for a < V — 2 for C"^~^'^ domains, let us note that the main ingredients of the results in 
Costabel (1988) are the near-optimal trace theorem for Lipschitz domains, which appears 
in Costabel (1988) and Ding (1996), and a certain regularity result for the Poincare-Steklov 
operator for Lipschitz domains, which appears, e.g., in McLean (2000). The relevant version 
of the trace theorem has been proved for C^^^'^ domains in Kim (2007), see also Marschall 
(1987). For the regularity of the Poincare-Steklov operator, the author has not been able 
to locate in the literature a result strong enough to give the boundedness for a < z^ — ^ , 
although there are results, e.g., in McLean (2000), that imply a = ly — 1 for standard 
boundary integral operators. 

The proof of Theorem 6.1 is divided into several lemmas that follow. The main idea 
is to decompose Tv into the part that is in Sp, which we call the low frequency part, and 
its complement, which we call the high frequency part. The low frequency part is readily 
handled by either the standard inverse estimates or the new inverse estimates from Dahmen, 
Faermann, Graham, Hackbusch, and Sauter (2004). To treat the high frequency part, we 
introduce a wavelet basis for the complement of Sp, and as naturally suggested by the 
techniques we use, the high frequency part is further decomposed into terms corresponding 
to far-field, near-field, and local interactions. Please be warned that these names are only 
suggestive in that, e.g., the local interaction terms may contain interactions between two 
wavelets with non-overlapping supports, although they cannot be too far apart. 

Our main analytic tool is a locally supported wavelet basis for the energy space H'^, 
with the dual multiresolution analysis based on piecewise polynomial- type spaces. More 
specifically, we assume that there is a Riesz basis ^ = {V'aIagV of H^ of wavelet type, 
whose dual, denoted by ^ = {V'aIagV) is focally supported piecewise polynomial wavelets, 
where V is a countable index set. Now we expand on what we mean exactly by the various 
adjectives such as "wavelet type" that characterize the bases ^ and ^. 

The collections ^ C H^{il) and ^ C H^^{fl) are biorthogonal: (V'AjV'/i) = ^x^J.^ ^-nd ^^^ 
Riesz bases for their corresponding spaces, meaning that 

\\T.XGS7Vxij\\\fjt(^n)^\\ivx)x\\i^\7), Ea6V^A^A ^_ ~||(ua)a||£2(v), (174) 

for any sequence {vx)x G £^(V). Here the notation X ^ Y means Y < X < Y. Each 
wavelet ipx or ipx has a scale, which is encoded by the function | • | : V — )• N. We say that 
if^x and Ipx have the scale 2~l ', which is justified by the locality properties 

diam(suppV'A) < 2~l^l, diam(supp^A) < S^'^L (175) 

We will also assume that the wavelet supports are locally finite, in the sense that 

#{A : |A| = e, B{x, 2-') n supp Va / 0} < 1, (176) 
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and similarly for the dual wavelets, where the bounds do not depend on ^ S N and x G il, 
and B{x,p) = {y G il : dist(a;,y) < p}. An immediate consequence of this property is that 
#{AGV:|A|=n<2"^. 

We assume that the wavelets have the so-called cancellation property of order d € N, 
saying that^ there exists a constant rj > 0, such that for any p G [l,oo], for all continuous, 
piecewise smooth functions v on Pq and A G V, 

\{v,ijjx)\ < 2"'^'^^"^^*^'^^max|t;|^^d,p(B(,upp^^^2-l^lr,)nr)> (177) 

where for ^ C M^ and e > 0, B{A, e) := {y £R^ : dist(A, y) < e}. 

Furthermore, we assume that for all r G [— d, 7), s < 7, necessarily with |s|, \r\ < v, 

\\w\\Hr[tt) < 2^^'^-'^\\w\\Hs(n), for w G spanj^-A : |A| = £}, (178) 

with 7 = sup{s : ^ C H^{Q)}, and similarly for the dual wavelets, with 7 and d replaced 
by 7 and d, respectively. 

We assume that the norm equivalence 

I^^-v-a^a||' ^-E2^^^^*^'"|-aP, (179) 

^ ' agv 

is valid for s G (0, d) D (0,7). 

For A G V, we define fi;^ = i7 n Bx, with Bx an open ball with diam(i?;^) < 2^1'^!, 
containing both suppV'A and supp^A- Thus Qx can be thought of as a common support of 
ipx and ipx- We then define 

A = {A : f7A n r / and diam(r2A) > Sh^ for some r G P}, (180) 

where (5 > is a constant so small that r^A n r / and A A imply fix C uj{t) for t G P, 
and assume that 

\\v\\H^in) ^Whp'^vW, f G S'a := spanlV-A : A G A}, (181) 

for s G [0, u]. By a duality argument (test v against w ~ /ip f ) this implies 

\\hpv\\ < \\v\\H-s(n), V G 5a, (182) 

for s G [0,1^]. 

Remark 6.4. Concrete examples of wavelet bases satisfying all our assumptions are given 
by the duals of the bases constructed in Dahmen and Schneider (1999). On triangulations 
over a Lipschitz polyhedral surface, one can also use the construction in Stevenson (2003). 
The only reason for insisting on polynomial dual wavelets is that in the proof of Lemma 
6.5 below, we use the inverse estimate (182), which is a consequence of (181). The latter 
estimate (hence both) can be proven by adapting the techniques from Dahmen et al. (2004). 

^ For notational convenience we interchanged the roles of d and d as compared to, e.g., Stevenson (2004). 
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Now that we have settled on our main tool, we can start with Proof of Theorem 6.1. 
Let V be as in the theorem. Then first we estimate the part of Tv that is in 5a. We define 
the projection operator Q\ : H~^{Q) — ;• 5a by 

Qa ^ wxijjx = ^ wx'ipx. (183) 

Aev AeA 

In what follows, we will abbreviate the Sobolev norms as || • ||s,t^ = || • ||Hs(tj) ^ind 

II ■ lis = II • llH»(f7)- 

Lemma 6.5 (Low frequency). Let s > 0, and suppose that either t>Oors + t>0. Then 
we have 

Y.hl(^-^'^\\Qj,Tv\\l^^,-^<\\vf, for vGH\n). (184) 

zeNp 

Proof. First let t < 0, and therefore s + 1 > 0. Then we have 

Y^ hl'^^^'^WQATvWl^^,^ < ^ IIQaT.;^.,^,) < WQaTvWI, < \\Tvf_, < \\v\\l (185) 

z£Np z£Np 

where we have used in succession a standard inverse estimate, the super-additivity of the 
Sobolev norms (3), the stability of Q\ in H~^, and the boundedness of T. 
For the case t > 0, a standard inverse estimate gives 

j; hl^^+'^QATvWi^^,^ < j; hfWQATvWl^,^ = Y. ||/.*pQaTH|^(,) 
zeNp zeNp zeNp (186) 

< Wh'pQATvf. 

At this point we employ the inverse estimate (182), to get 

\\h'pQATv\\ < WQaTvU < \\TvU < \\v\\t, (187) 

concluding the proof. □ 

What remains now is to bound (/ — Qa)Tv, which consists of only high frequency 
wavelets compared to what is in QaTv G Sa- To this end, for A G A^ := V \ A, let us 
define £x by 

2"^^ = max{hr : r G P, r n S^a / 0}, (188) 

so that in light of the norm equivalence (179), we can write 

y: hi^^^mi-QA)Tv\\i^^^^ < y: hi^'^'^ E 22i^i(^+*)|(r.),|2 

z&Np zeNp {AeA=:riAnaj(z)7^0} 

< ^ E 2"2^^(^+*)22|A|(^+*)|(r.),p (189) 



zeNp {X(^A'=:Qxr\io{z)j^0} 
AeA= 



^ y^ 22(|A|-^A)(s+t)|(2^^) 
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where {Tv)\ = {Tv, tp\) is the coordinate of Tv with respect to ipx^ and in the last step we 
have taken into account the fact that each Q\ intersects with only a uniformly bounded 
number of stars u){z). Hence our aim is to bound the last expression in (189) by ||(fA)A||^2 
for V G Sp, where v\ = {v,tpx)- 

Before dealing with nonlocality of T, let us focus on the local properties. To this end, 
with Sp^z = {v € Sp : V = outside uj^{z)}, define Qz to be the L^-orthogonal projec- 
tor onto Sp^z in case of piecewise constants, and otherwise to be the quasi-interpolation 
operator onto Sp^z, as discussed around (21). Here uj^{z) = uj{uj^~^{z)) for k > 2. 

Lemma 6.6 (Local interactions). Let s > and let t + s < a. Then we have 

Y,hf^+'^\\{I-QA)TQzv\\l^^^^<\\vf, for v e Sp. (190) 

z£Np 

Proof. Replacing v with QzV in (189), we get 

E hl^'^'^\\iI-QA)TQ.v\\l^iz) < E 2^(l^l-^^)(^^*)|(rQ..)A|^ (191) 

zeNp AeA<: 

In order to bound this, we view it as the (squared) £^(A'^)-norm of a sequence, and we 
test the sequence against a general sequence {wx)x G ^^(A'^). As a preparation to this, for 
w G spanjf/'A : |A| = ^}, we have 

\{w,TQzv)\ < \\w\\t-jTQzV\\^t+a < 2-^^\\w\\t\\Qzv\\t+a 

where £z is defined by 2"^^ = hz- Then assuming that t < 0, for w = ^j^tt'AV'A with 
{w\)\ G ^^(A'^), we infer 






E wxip\,TQzV 



^ E E 2(^-^^)(^+*-)+^^* E i^aP I ii^^iL3(.) 
^ E 2'^* E i^^aI'I ii^l^w ^ ii^iui^p ^11 ^ ii^iitii^iit' (193) 

zeNp \{\<^A'':nxr\u{z)^0} J 

establishing the lemma. 

The case i > is proven similarly, by using 

\{w,TQzv)\ < 2-'^'+'^'^\\w\\t\\v\U^^^z), (194) 

instead of (192). □ 
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To deal with nonlocality of T, we will employ the Schur test, whose proof we include 
here for convenience. 

Lemma 6.7 (Schur test). Let A and M be countable sets, and let £p{A,w) be the weighted 
iP -space with the norm 

\Mep{A,w) = \\{vxwx}x(^A\\ep(A), (195) 

where {w\}x^\ is a given positive sequence of weights. For a matrix with entries T\^, 
A G A, ^ G M, we consider its boundedness as a linear operator T : l!P{M) — t- I!P{K,w). 
Suppose that there exist two positive sequences {9\\\^\ and {u:^}^^m, cmd two numbers 
a, /3 G M such that 

/leM AeA 

where 1 < p < oo and q satisfies - -\ — = 1. Then we have 

\\T\\£P(M)-^ep{A,w) < "/3- (197) 

Proof. Let v G P'{M). Then we have 

1 1 

<3 / \ P 






pi^H'Kr 



(198) 






Mi^i^^r 



and so 



.mKI^mI" 



r-ii?.(AH = E i(^-)aIM < «^E^>A E i^^, 

AeA AeA fieM 

= aPY ^i^/^i" E i^A,i^>^ < a^ E ^;:'\v,\'P'ujp (199) 

/^eAf AeA /xeA/ 

— nP8P\\v\\P 

estabhshing the boundedness of T : F(M) -^ ^^(A, w). □ 

With the Riesz bases ^ C H\Q) and § C H^\Q) at hand, the operator T : H^Q) -^ 
H^^{Q) can be thought of as the bi-infinite matrix T : £^(V) — t- ^^(V) with the elements 
^A^ = (V'A)2^'0/i)- The following estimate on these elements, established in Stevenson 
(2004), will be crucial: 

|rA.| = K^A,T^.)|<^(3^;^ (200) 
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where 

5{X, fi) = 2™'^^l^l'l^l>dist(f]A, %)■ (201) 

See also Schneider (1998) and Dahmen and Stevenson (1999) for earher derivations for less 
general cases. 

The following result shows that the far-field terms behave rather well. Note that the 
condition d -|- t > is immaterial to us since we can choose d at will. 

Lemma 6.8 (Far-field interactions). Define the operator F : i/*(il) — t- H^''(Q) with matrix 
elements Fx^ = T\^ for A G A'^ and ^u G V satisfying d\si{Q.\,^^) > (5max{2~^^, 2~l^l}, 
and Fxfj^ = otherwise. Let s -|- t > and d + t > 0. Then we have 

^ 22(l^l-^A)(s+t)|(^^)^|2< ||y||2 j^^ veH\Q). (202) 

Proof We apply the Schur test (Lemma 6.7) with 0l = 2-l-^l"/2-{|A|-£A){rf+t)^ ^;2 ^ 2-l'^l"/2, 
and wx = 2(l''*l^^^)('^+*). First we shall bound 

(-IIAI- 1/2 \ "+2*+2d 

by a multiple of Q\. In the inner sum, we must sum over all [i such that (5(A,//) > 2l^l^^^ 
when \\x\=m> |A|, and such that d{\,fj,) > 2™'*^i0'''"~^^i when \n\ = m < |A|. By locality 
of the wavelets and the Lipschitz property fi, for A G V, m € N and /? > 0, we have 

J2 ^(^^^)-(n+/3) <^~/32"'^a^{0.™-|A|}^ (204) 

{li:\fi\=m,5{\,ii)>R} 

which appears, e.g., in Stevenson (2004). This gives 

V j(-;^^^>)-(n+2t+2d) < 2-(|A|-£A)(2rf+2t)+(m-|A|)n/2^ (205) 

{/^:|At|=m,<5(A,/^)>2l^|-*A} 

for m> |A|, and 

V 5(A,/i)"("+^*+^"') < 2"(l^l"^^)('^+*)+('""l^l)'^/^ (206) 

{/i:|/i|=m,<5(A,At)>2™="'-tO'™-^A}} 

for m < |A|. By using the preceding estimates, we conclude 

V |Fa |2"I^I"/2 < y^ 2-(i-^i-^>')('^+*)-(i^i-")('^+*)-i^i"/2 

AieV m<|A| 

+ y^ 2-{\M-ex){2d+2t)-{m-\X\){d+t)-\X\n/2 (207) 

m>|A| 
< 2-|A|n/2-(|A|-£;,)(d+t)_ 
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Now we shall bound 

by a multiple of uJj^. By construction, for A G A'^ and /x G V with Fx^ / 0, we have 
2-mm{£^,|^|} < (iist(f]A, ^^l), implying that 

2-^^ <dist(f7A,fi;.), and dist(J]A,^M)^2-^"^'^^l^l'l^l>. (209) 

In particular, the second estimate tells us that 6{X,fj.) > 1 in the (inner) sum in (208). So 
for I < |/i|, we estimate the inner sum as 

\\\=i \\\=£ 

/ n \ e\ /n\ n+2t+2d 



{\x\=e,s(x,f,)>i} \ "■ '^' 

< 2-in/22-i\i^\-e){n/2+t+d) ^ 2-(l/^l-^)(*+'=')2~l'^l"/ 



2 



where we have used i > £\ and d + t > in the first inequality, and s + i > in the second. 
Assuming that 2s + t — d < 0, for i > \fj,\, we estimate 



(211) 



Y^ 2{e~e,)i2s+t~d)2-en/2^j,^^^ < ^ 2-^"/2|Fa/.| 

|A|=£ |A|=^ 

< 2-in/22-{i-\^l\){n/2+t+d)2ii-\^i\)n 

= 2-{MfJ'\)(d+t)2-~Mn/2^ 

Now if 2s + i — d > 0, for i > \fi\, we have 

y^ 2(^-^^)(2'^+*-'^)2"^"/2|Fa^| < V 2('^~I''I)(2'+*"'^)(5(A,^)2'+*"''2"^"/2|Fa^| 

|A|=£ |A|=£ 

< 2{^-|/^l)(2s+i-'^)2-M22-{^-|Ml){n/2+t+d)2(^-|Ml)'^ 
= 2-(^-|H)(2rf-2s)2-|M|n/2^ 

From the geometric decay of the preceding estimates, we infer 

y^ 2(I^I~^^)(2«+*"'^)2-I^I"/2|Fa I < 2^1^!"/^ (213) 

AeA'^ 

which completes the proof. □ 



(212) 
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For the remaining terms, i.e., for the near-field terms, we employ the simple estimate 

{i'X.TlP^) < Uy^Wt-ATi^^W-t+a < Ux\\t-Ai'A\t+. < 2-|^l-+H-. (214) 

We will see that this estimate gives sub-optimal results, that in general require the manifold 
fi to be smoother than Lipschitz. There exist sharper estimates in the literature, cf. 
Stevenson (2004); Dahmen et al. (2006), that exploit the piecewise smooth nature of the 
wavelets. However, the author was not able to make use of them to get better results. 
The best attempts so far by using the estimates from Dahmen et al. (2006) resulted in 
logarithmic divergences, and the estimates from Stevenson (2004) in general need Vl to be 
smoother than Lipschitz, thus do not seem to give improvements in this regard. 

Lemma 6.9 (Near-field interactions). Define the operator N : H*{Q) — )• H~^{Q,) with 
matrix elements Nx^ = T\^ for A G A^ and /i € V satisfying disi{Vlx,Q.y) < 2~l'^l and 
\lA ^ ^A? CLnd N\^ = otherwise. Let 2a > n. Then we have 

^22(I^M^)(^+t)|(Ar„),|2<||,||2 f^^ ^^sp. (215) 

AeA= 

Proof We apply the Schur test (Lemma 6.7) with Sf = 2^^^-''l^l, w^ = 2^e-a)M^ ^^^ 
wx = 2(l^l-'^^)(^+*), where e = 2a - n > 0. We have 

Y^ |iVAM|2^^"''^'''' < Yl 2"'"l^l+^" < 2^^^-''l^l, (216) 

and 

y^wl6l\Nx^\< V \:^22(^~^")(''+*)2=^""''^2"''(^"l''l)2"(^"^") 

AeA= {reP:dist(r,n,i)<2-|Ml}£>£r 

< y^ 2^^^~'^^^2"'^(^^"i^i) 

{reP:dist(r,n^)<2-lMl} (217) 

= Y 2'"l''l-"^- < Y 2'"l''lvol(r) 

{reP:dist(T,n,i)<2-|Ml} {reP:dist(T,Cp)<2-|Ml } 

< 2('^^")l^l = 2(^~°")l'^l 
which establishes the proof. □ 

The proof of Theorem 6.1 is complete, as choosing 6 > small enough, the combination 
of the local, far-field and near-field interactions accounts for the high frequency part (/ — 
Q\)Tv, while the low frequency part QaTv is taken care of by Lemma 6.5. 
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7 Concluding remarks 

In this paper, we proved geometric error reduction for three kinds of adaptive boundary 
element methods for positive, negative, as well as zero order operator equations, without 
relying on a saturation- type assumption. In fact, several types of saturation assumptions 
follow from our work as a corollary. Moreover, bounds on the convergence rates are obtained 
that are in a certain sense optimal. 

We established several new global- and local discrete bounds for a number of residual- 
type error estimators for positive, negative, as well zero order boundary integral equations, 
including the estimators from Carstensen, Maischak, and Stephan (2001), Carstensen, 
Maischak, Praetorius, and Stephan (2004), and Faermann (2000, 2002). Some of our 
bounds contain oscillation terms, that give useful estimates on how far the current mesh 
is from saturation. In order to handle the oscillation terms, which turned out to be much 
more involved than the corresponding affairs in the finite element context, we introduced an 
inverse-type inequality involving boundary integral operators and locally refined meshes. 
Our proof of the inequality in general requires the underlying surface T to be C^'^ or 
smoother, but for open surfaces it allows the boundary of T to be Lipschitz. So in general, 
polyhedral surfaces are ruled out, which is likely an artifact of the proof, and will hopefully 
be rectified in the near future. 

The current work gives rise to its fair share of open problems, and re-emphasizes some 
existing ones. The following is an attempt at identifying the most pressing of them. 

• to prove the inverse-type inequality in §6 for polyhedral surfaces. 

• to characterize the approximation classes associated to the proposed adaptive BEMs. 

• to extend the analysis to transmission problems, and adaptive FEM-BEM coupling. 

• complexity analysis, i.e., the problem of quadrature and linear algebra solvers. 

• convergence rate for adaptive BEMs based on non-residual type error estimators. 
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